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Recent experiments on triangular lattice organic Mott insulators have found evidence for a 2D spin liquid 
in close proximity to the metal-insulator transition. A Gutzwiller wavefunction study of the triangular lattice 
Heisenberg model with a four-spin ring exchange term appropriate in this regime has found that the projected 
spinon Fermi sea state has a low variational energy. This wavefunction, together with a slave particle gauge 
theory analysis, suggests that this putative spin liquid possesses spin correlations that are singular along surfaces 
in momentum space, i.e. "Bose surfaces". Signatures of this state, which we will refer to as a "Spin Bose-Metal" 
(SBM), are expected to be manifest in quasi-lD ladder systems: The discrete transverse momenta cut through 
the 2D Bose surface leading to a distinct pattern of ID gapless modes. Here, we search for a quasi- ID descendant 
of the triangular lattice SBM state by exploring the Heisenberg plus ring model on a two-leg triangular strip 
(zigzag chain). Using DMRG supplemented by variational wavefunctions and a Bosonization analysis, we map 
out the full phase diagram. In the absence of ring exchange the model is equivalent to the Ji — J2 Heisenberg 
chain, and we find the expected Bethe-chain and dimerized phases. Remarkably, moderate ring exchange reveals 
a new gapless phase over a large swath of the phase diagram. Spin and dimer correlations possess singular 
wavevectors at particular "Bose points" (remnants of the 2D Bose surface) and allow us to identify this phase 
as the hoped for quasi-lD descendant of the triangular lattice SBM state. We use Bosonization to derive a 
low energy effective theory for the zigzag Spin Bose-Metal and find three gapless modes and one Luttinger 
parameter controlling all power law correlations. Potential instabilities out of the zigzag SBM give rise to other 
interesting phases such as a period-3 Valence Bond Solid or a period-4 Chirality order, which we discover in 
the DMRG. Another interesting instability is into a Spin Bose-Metal phase with partial ferromagnetism (spin 
polarization of one spinon band), which we also find numerically using the DMRG. 



I. INTRODUCTION 

A promising regime to search for elusive 2D spin Uquids 
is in the proximity of the Mott metal-insulator transition. In 
such "weak Mott insulators" significant local charge fluctua- 
tions induce multi-spin ring exchange processes which tend to 
suppress magnetic or other types of ordering. Indeed, recent 
experiments''^ on the triangular lattice based organic Mott in- 
sulator K-(ET)2Cu2(CN)3 reveal no indication of magnetic 
order or other symmetry breaking down to temperature sev- 
eral orders of magnitude smaller than the characteristic ex- 
change interaction energy J w 250 K. Under pressure the k- 
(ET)2Cu2(CN)3 undergoes a weak first order transition into a 
metallic state, while at ambient pressure it has a small charge 
gap of 200 K, as expected in a weak Mott insulator. Ther- 
modynamic, transport, and spectroscopic experiments''^'"* all 
point to the presence of a plethora of low energy excitations 
in the k-(ET)2Cu2(CN)3, indicative of a gapless spin liquid 
phase. Several authors have proposed^"^ that this putative 
spin liquid can be described in terms of a Gutzwiller-projected 
Fermi sea of spinons. 

Quantum chemistry calculations suggest that a one-band 
triangular lattice Hubbard model at half filling is an appropri- 
ate theoretical starting point to describe k-(ET)2Cu2(CN)3.''^ 
Variational studies of the triangular lattice Hubbard model^ 
find indications of a non-magnetic spin liquid phase just on the 
insulating side of the Mott transition. Moreover, exact diag- 



onalization studies of the triangular lattice Heisenberg model 
show that the presence of a four-site ring exchange term ap- 
propriate near the Mott transition can readily destroy the 120° 
antiferromagnetic order. One of us^ performed variational 
wavefunction studies on this spin model and found that the 
Gutzwiller-projected Fermi sea state* has the lowest energy 
for sufficiently strong four-site ring exchange interactions ap- 
propriate for the K-(ET)2Cu2(CN)3. 

Despite these encouraging hints, the theoretical evidence 
for a spin liquid phase in the triangular lattice Hubbard model 
or Heisenberg spin model with ring exchanges is at best sug- 
gestive. Variational studies are biased by the choice of wave- 
functions and can be notoriously misleading. Exact diagonal- 
ization studies are restricted to very small sizes, which is es- 
pecially problematic for gapless spin liquids. Quantum Monte 
Carlo fails due to the sign problem. The density matrix renor- 
malization group (DMRG) can reach the ground state of large 
ID systems, but capturing the highly entangled and non-local 
character of a 2D gapless spin liquid state is a formidable chal- 
lenge. Thus, with new candidate spin liquid materials and 
increasingly refined experiments available, the gap between 
theory and experiment becomes ever more dire. 

Effective field theory approaches such as slave particle 
gauge theories or vortex dualities, while unable to solve any 
particular Hamiltonian, do indicate the possibility of stable 
gapless 2D spin liquid phases. Such gapless 2D spin liquids 
generically exhibit spin correlations that decay as a power 



law in space, perhaps with anomalous exponents, and which 
can oscillate at particular wavevectors. The location of these 
dominant singularities in momentum space provides a con- 
venient characterization of gapless spin liquids. In the "al- 
gebraic" or "critical" spin liquids""''^ these wavevectors are 
limited to a finite discrete set, often at high symmetry points 
in the Brilloin zone, and their effective field theories can of- 
ten exhibit a relativistic structure. But the singularities can 
also occur along surfaces in momentum space, as they do in 
the Gutzwiller-projected spinon Fermi sea state, the 2D Spin 
Bose-Metal (SBM) phase. It must be stressed that it is the 
spin (i.e., bosonic) correlation functions that possess such sin- 
gular surfaces - there are no fermions in the system - and the 
low energy excitations cannot be described in terms of weakly 
interacting quasiparticles. It has been proposed recently'^ 
that a 2D "Boson-ring" model describing itinerant hard core 
bosons hopping on a square lattice with a frustrating four- 
site term can have an analogous liquid ground state which we 
called a d-wave Bose liquid (DBL). The DBL is also a Bose- 
Metal phase, possessing a singular Bose surface in momentum 
space. 

Recently we have suggested'* '^ that it should be possi- 
ble to access such Bose-Metals by systematically approach- 
ing 2D from a sequence of quasi-lD ladder models. On a 
ladder the quantized transverse momenta cut through the 2D 
surface, leading to a quasi- ID descendant state with a set of 
low-energy modes whose number grows with the number of 
legs and whose momenta are inherited from the 2D Bose sur- 
faces. These quasi-lD descendant states can be accessed in 
a controlled fashion by analyzing the ID ladder models us- 
ing numerical and analytical approaches. These multi-mode 
quasi- ID liquids constitute a new and previously unantici- 
pated class of quantum states interesting in their own right. 
But more importantly they carry a distinctive quasi- ID "fin- 
gerprint" of the parent 2D state. 

The power of this approach was demonstrated in Ref. 16 
where we studied the new Boson-ring model on a two-leg lad- 
der and mapped out the full phase diagram using the DMRG 
and ED, supported by variational wavefunction and gauge the- 
ory analysis. Remarkably, even for a ladder with only two 
legs, we found compelling evidence for the quasi- ID descen- 
dant of the 2D DBL phase. This new quasi- ID quantum state 
possessed all of the expected signatures reflecting the parent 
2D Bose surface. 

In this paper we turn our attention to the triangular lat- 
tice Heisenberg model with ring exchange appropriate for 
the K-(ET)2Cu2(CN)3 material. In hopes of detecting the 
quasi- ID descendant of the triangular lattice Spin Bose-Metal 
(Gutzwiller-projected spinon Fermi sea state), we place this 
model on a triangular strip with only two legs shown in Fig. 1 . 
The all-important ring exchange term acts around four-site 
plackets as illustrated; we also allow different Heisenberg ex- 
change couplings along and transverse to the ladder. 

It is convenient to view the two-leg strip as a Ji — J2 chain 
(studied extensively before'^ '^) with additional four-spin ex- 
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FIG. 1: Top: Heisenberg plus ring exchange model on a 2-leg tri- 
angular strip. Bottom: Convenient representation of the model as 
a Ji — J2 chain with additional four-site terms; the Hamiltonian is 
written out in Eq. (1). 



changes. The Hamiltonian reads 

H = ^ [ 2Ji5(x) • S{x + 1) + 2J2Six) ■ S{x + 2) 

X 

+ H.C.)]. (1) 

The four-spin operators act as P1234 : |cri, (T2, 0-3, (T4) 
|(T4, (Ji, (72, (T3), -P/234 — (^1234)^^- We attack this model 
using a combination of numerical and analytical techniques - 
DMRG, exact diagonalization (ED), variational Monte Carlo 
(VMC), as well as employing bosonization to obtain a low en- 
ergy effective field theory from the slave particle-gauge for- 
mulation (and/or from an interacting electron Hubbard-type 
model). Our key findings are summarized in Fig. 2 which 
shows the phase diagram for antiferromagnetic couplings Ji, 
J2, and luring- For K,:iag — > 0, the Ji — J2 model has the 
familiar ID Bethe-chain phase for J2 < 0.24 Ji and period- 
2 Valence Bond Solid (VBS-2) for larger J2. For A'ring ^ 
0.2 Ji, new physics opens up. In fact, Klironomos et al?^ con- 
sidered such Ji — J2 — A'ling model motivated by the study 
of Wigner crystals in a quantum wire.^' Using ED of systems 
up to L = 24, they found an unusual phase in this intermedi- 
ate regime, called "4P" in Fig. 8 of Ref. 20, but it had proven 
difficult to clarify its nature. We identify this region as a de- 
scendant of the triangular lattice Spin Bose-Metal phase (or 
further derivatives of the descendant as discussed below). 

A caricature of the zigzag Spin Bose-Metal is provided 
by considering a Gutzwiller trial wavefunction construction 
on the two-leg strip. The 2D SBM is obtained by letting 
spinons hop on the triangular lattice with no fluxes and then 
Gutzwiller-projecting to get a trial spin wavefunction. So here 
we also take spinons hopping on the ladder with no fluxes, 
which is ti — 1-2 hopping in the ID chain language that we 
mainly use. For t2 < 0.5ii, the mean field state has one 
Fermi sea segment spanning [— 7r/2, 7r/2] (spinons are at half- 
filling), and the Gutzwiller projection of this is known to be 
an excellent state for the Bethe chain. On the other hand, for 
t2 > 0.5ti, the spinon band has two Fermi seas as shown in 
Fig. 3. The Gutzwiller projection of this is a new phase that 
we identify as a quasi- ID descendant of the triangular lattice 
Spin Bose-Metal. The wavefunction has one variational pa- 
rameter 12/ ti, or, equivalently, the ratio of the two Fermi sea 
volumes. Using this restricted family of states, our VMC ener- 
getics study of the Ji — J2 — ATring model finds three regimes 
broadly delineated by solid lines in Fig. 3 for larger /Ving: 
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FIG. 2: (Color online) Phase diagram of the ring model Eq. (1) deter- 
mined in the DMRG using system sizes L = 48 — 96. Filled squares 
(red) denote Bethe-chain phase. Open squares (with black outlines) 
denote Valence Bond Solid with period 2. Open circles (blue) denote 
Spin Bose-Metal. Open circles with crosses denote where the DMRG 
has difficulties converging to singlet for the larger sizes, but where 
we still think this is spin-singlet SBM. Star symbols denote points 
where the ground state appears to have true non-zero spin (for all 
points here, the magnetization is smaller than full polarization of the 
smaller Fermi sea in the SBM interpretation). Filled diamonds (ma- 
genta) denote VBS with period 3. Our identifications are ambiguous 
in the lower VBS-3 region approaching VBS-2. Lines indicate phase 
boundaries determined in VMC using spin-singlet wavefunctions de- 
scribed in the text (we also used appropriate dimerized wavefunc- 
tions for the VBS states). A detailed study of a cut ivring/ Ji = 1 is 
presented in Sec. Ill (cf. Fig. 8). 



i) In the Bethe-chain regime the optimal state has one Fermi 
sea. ii) For sufficiently large A'ling and upon increasing J2, 
we enter a different regime where it is advantageous to start 
populating the second Fermi sea. As we further increase J2 
moving away from the Bethe-chain phase, we gradually trans- 
fer more spinons from the first to the second Fermi sea. This 
whole region is the SBM. Hi) Finally, at still larger J2, the vol- 
umes of the two Fermi seas become equal, which corresponds 
to t2/ti 00, i.e., decoupled-legs limit. 

The DMRG is the crucial tool that allows us to answer 
how much of this trial state picture actually holds in the 
Ji — J2 — Kj-ing model. Fig. 2 shows all points that were 
studied using the DMRG and their tentative phase identifica- 
tions by looking at various ground state properties. Remark- 
ably, in a broad-brush sense, the three regimes found in VMC 
for Kring > 0.2 (one spinon Fermi sea, two generic Fermi 
seas, and decoupled legs) match quite closely different qual- 
itative behaviors found by the DMRG study and marked as 
Bethe-chain, SBM, and VBS-2 regions. Here we note that the 
decoupled-legs Gutzwiller wavefunction is gapless and does 
not have VBS-2 order, but it is likely unstable towards opening 
a spin gap;"^''** still, it is a good initial description for large J2- 
On the other hand, away from the decoupled-legs limit, we ex- 
pect a stable gapless SBM phase. The DMRG measures spin 
and dimer correlations, and we identify the SBM by observing 
singularities at characteristic wavevectors that evolve continu- 
ously as we move through this phase - these are the quasi- ID 
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FIG. 3: Spinon dispersion for t2 > 0.5ti showing two occupied 
Fermi sea segments (here and throughout we use the ID chain lan- 
guage, see bottom Fig. 1). Gutzwiller projection of this is the zigzag 
SBM state at the focus of this work. 



"Bose points" (remnants of a 2D Bose surface). The singular 
wavevectors are reproduced well by the VMC, although the 
Gutzwiller wavefunctions apparently cannot capture the am- 
plitudes and power law exponents. 

An effective low energy field theory for the zigzag SBM 
phase can be obtained by employing Bosonization to ana- 
lyze either a spinon gauge theory formulation or an inter- 
acting model of electrons hopping on the zigzag chain. In 
the latter case we identify an Umklapp term which drives the 
two-band metal of interacting electrons through a Mott metal- 
insulator transition. The low energy Bosonized description of 
the Mott insulating state thereby obtained is identical to that 
obtained from the zigzag spinon gauge theory. In the interact- 
ing electron case, there are physical electrons that exist above 
the charge gap. On the other hand, in the gauge theory the 
"spinons" are unphysical and linearly confined. 

The low energy fixed point theory for the zigzag Spin Bose- 
Metal phase consists of three gapless (free Boson) modes, 
two in the spin sector and one in the singlet sector (the lat- 
ter we identify with spin chirality fluctuations). Because of 
the SU(2) spin invariance, there is only one Luttinger param- 
eter in the theory, and we can characterize all power laws 
using this single parameter. The dominant correlations oc- 
cur at wavevectors 2kpi and 2kF2 connecting opposite Fermi 
points, and the power law can vary between x^^^^ and x^^ 
depending on the value of the Luttinger parameter We un- 
derstand well the stability of this phase. We also understand 
why the Gutzwiller-projected wavefunctions, while capturing 
the singular wavevectors, are not fully adequate - our trial 
wavefunctions appear to be described by a specific value of 
the Luttinger parameter that gives x~^/^ power law at 2fci?i 
and 2kF2- The difference between the DMRG and VMC in 
the SBM phase is qualitatively captured by the low energy 
Bosonized theory. 

The full DMRG phase diagram findings are in fact much 
richer. Prominently present in Fig. 2 is a new phase occur- 
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ring inside the SBM and labeled VBS-3. This has period-3 
Valence Bond Solid order "dimerizing" every third bond and 
also has coexisting effective Bethe-chain-like state formed by 
non-dimerized spins (see Sec. V A and Fig. 14 for more ex- 
planations). A careful look at the SBM theory reveals that at 
a special commensuration where the volume of the first Fermi 
sea is twice as large as that of the second Fermi sea, the SBM 
phase can be unstable gapping out the first Fermi sea and pro- 
ducing such VBS-3 state. 

Another observation in Fig. 2 is the possibility of develop- 
ing a partial ferromagnetic (FM) moment in the SBM region 
labeled "partial FM" to the left of the VBS-3 phase. We do 
not understand all details in this region. In the SBM further 
to the left, we think the ground state is spin-singlet, which is 
what we find from the DMRG for smaller system sizes up to 
L = 48. However, the DMRG already has difficulties con- 
verging to the spin-singlet state for the larger system sizes 
L ~ 96. In the partial FM region, it seems that the ground 
state has a small magnetization. Given our SBM picture, it is 
conceivable that one or both spinon Fermi seas could develop 
some spin polarization. The most likely scenario is for the 
polarization to first appear in the smaller Fermi sea since it is 
more narrow in energy (more flat-band-like). The total spin 
that we measure in the partial FM region in Fig. 2 is smaller 
than what would be expected from a full polarization of the 
second Fermi sea, and it is hard for us to analyze such states. 

To check our intuition, we have also considered a modified 
model with an additional third-neighbor coupling J3 which 
can be either antiferromagnetic or ferromagnetic tailored to 
either suppress the ferromagnetic tendencies or to reveal them 
more fully. We have studied this model at ifring = Ji = 1, 
varying J2. With antiferromagnetic J3 = 0.5, we have in- 
deed increased the stability of spin-singlet states in the region 
to the left of the VBS-3. Interestingly, this study, which is 
not polluted by the small moment difficulties, also revealed a 
new spin-gapped phase near the VBS-3. The new phase has a 
particular period-4 order in the spin chirality, and we can un- 
derstand the occurrence as an instability at another commen- 
suration point hit by the singular wavevectors as they vary in 
the SBM phase (see Sec. IV E for details). Turning now to 
ferromagnetic J3 = —0.5, we have found a more clear ex- 
ample of the partial ferromagnetism where the ground state 
is well described by Gutzwiller-projecting a state with a fully 
polarized second Fermi sea and an unpolarized first Fermi sea. 

The paper is organized as follows. To set the stage, in 
Sec. II we develop general theory of the zigzag SBM phase. 
In Sec. Ill we present the DMRG study of the ring model that 
leads to the phase diagram Fig. 2. We consider carefully the 
cut at A'ring = Ji and provide detailed characterization of 
the new SBM phase. In Sec. IV we study analytically the 
stability of the SBM. We also consider possible phases that 
can arise as some instabilities of the SBM. This is done in 
particular to address the DMRG findings of the VBS-3 and 
Chirality-4 states, which we present in Sec. V. To clarify the 
regime to the left of the VBS-3 where the DMRG runs into 
convergence difficulties or small moment development, we 
also perturb the model with antiferromagnetic (Sec. VB) or 
ferromagnetic (Sec. VI) third-neighbor interaction J3 and dis- 



cuss partially polarized SBM. Finally, in Sec. VII we briefly 
summarize and suggest some future directions one might ex- 
plore. In Appendix A, within our effective field theory anal- 
ysis, we summarize the bosonization expressions for physical 
observables that are measured in the DMRG and VMC. In Ap- 
pendix B we provide details of the Gutzwiller wavefunctions 
that are used throughout in the VMC analysis. In Appendix C 
we summarize the DMRG results for the conventional Bethe- 
chain and VBS-2 phases. 

II. SPIN BOSE-METAL THEORY ON THE ZIGZAG STRIP 

Since a wavefunction does not constitute a theory, and can 
at best capture a caricature of the putative SBM phase, it is 
highly desirable to have a field-theoretic approach. The goal 
here is to obtain an effective low energy theory for the SBM 
on the zigzag chain. In 2D the usual approach is to decompose 
the spin operators in terms of an SU(2) spinor - the fermionic 
spinons: 

S = i/taa/3//3 ; flu - 1 . (2) 

In the mean field one assumes that the spinons do not inter- 
act with one another and are hopping freely on the 2D lat- 
tice. For the present problem the mean field Hamiltonian 
would have the spinons hopping in zero magnetic field, and 
the ground state would correspond to filling up a spinon Fermi 
sea. In doing this one has artificially enlarged the Hilbert 
space, since the spinon hopping Hamiltonian allows for un- 
occupied and doubly-occupied sites, which have no meaning 
in terms of the spin model of interest. It is thus necessary 
to project back down into the physical Hilbert space for the 
spin model, restricting the spinons to single occupancy. If 
one is only interested in constructing a variational wavefunc- 
tion, this can be readily achieved by the Gutzwiller projection, 
where one simply drops all terms in the wavefunction with 
unoccupied or doubly-occupied sites. The alternate approach 
to implement the single occupancy constraint is by introduc- 
ing a gauge field, a U{1) gauge field in this instance, that 
is minimally coupled to the spinons in the hopping Hamil- 
tonian. This then becomes an intrinsically strongly-coupled 
lattice gauge field theory. To proceed, it is necessary to resort 
to an approximation by assuming that the gauge field fluctua- 
tions are (in some sense) weak. In 2D one then analyzes the 
problem of a Fermi sea of spinons coupled to a weakly fluc- 
tuating gauge field. This problem has a long history, '^'^^"^^ 
but all the authors have chosen to sum the same class of dia- 
grams. Within this (uncontrolled) approximation one can then 
compute physical spin correlation functions, which are gauge 
invariant. It is unclear, however, whether this is theoretically 
legitimate, and even less clear whether or not the spin liq- 
uid phase thereby constructed captures correctly the universal 
properties of a physical spin liquid that can (or does) occur for 
some spin Hamiltonian. 

Fortunately, on the zigzag chain we are in much better 
shape. Here it is possible to employ Bosonization to analyze 
the quasi- ID gauge theory, as we detail below. While this still 
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does not give an exact solution for the ground state of any 
spin Hamiltonian, with regard to capturing universal low en- 
ergy properties it is controlled. As we will see, the low energy 
effective theory for the SBM phase is a Gaussian field theory, 
and perturbations about this can be analyzed in a systematic 
fashion to check for stability of the SBM and possible insta- 
bilities into other phases. 

As we will also briefly show, the low energy effective the- 
ory for the SBM can be obtained just as readily by starting 
with a model of interacting electrons hopping on the zigzag 
chain, i.e. a Hubbard-type Hamiltonian. If one starts with in- 
teracting electrons, it is (in principle) possible to construct the 
gapped electron excitations in the SBM Mott insulator Within 
the gauge theory approach, the analogous gapped spinon ex- 
citations are unphysical, being confined together with a linear 
potential. Moreover, within the electron formulation one can 
access the metallic phase, and also the Mott transition to the 
SBM insulator 



A. SBM via Bosonization of gauge theory 

We first start by using Bosonization^"""*- to analyze the 
gauge theory.^-'"-'^ Motivated by the 2D triangular lattice with 
ring exchanges, we assume a mean field state in which the 
spinons are hopping in zero flux. Here the spinons are hopping 
on the zigzag strip with near-neighbor and second-neighbor 
hopping strengths denoted ti and t2- This is equivalent to 
a strictly ID chain with likewise first- and second-neighbor 
hopping. The dispersion is 



^(fc) = -2ti cos(fc) - 2t2 cos(2A:) - fi . 



(3) 



For t2 > 0.5ti there are two sets of Fermi crossings at wave 
vectors ztkpi and ±A;f2 as shown in Fig. 3. Our conven- 
tion is that fermions near kpi and are moving to the 
right; the corresponding group velocities are vi, V2 > 0. The 
spinons are at half-filling, which implies kpi + kp2 ^ — 7r/2 
mod 2tt. 

The spinon operators are expanded in terms of continuum 
fields. 



a,P 



(4) 



with a = 1,2 denoting the two Fermi seas, a =t, j denoting 
the spin, and P ^ R/ L = ± denoting the right and left mov- 
ing fermions. We now use Bosonization,^""^- re-expressing 
these low energy spinon operators with Bosonic fields. 



with canonically conjugate boson fields: 



(5) 



[Paa{x),ipbl3{x')\ ^ [9aa{x),9bf3ix')] = , (6) 
[ipaaix),dbi3{x')] = inSabSaP Q{x - x') , (7) 

where Q{x) is the Heaviside step function. Here, we have in- 
troduced Klein factors, the Majorana fermions {ijaa, Vb/s} = 



25abSap, which assure that the spinon fields with different 
flavors anti-commute with one another. The (slowly vary- 
ing) fermionic densities are simply fj^^^fpaa = dx{Pfaa + 

eaa)/(2^). 

A faithful formulation of the physical system in this slave 
particle approach Eq. (2) is a compact U(l) lattice gauge the- 
ory. In 1h-1D continuum theory, we work in the gauge elimi- 
nating spatial components of the gauge field. The imaginary- 
time bosonized Lagrangian density is then: 



-E 

2tt ^ 



(8) 



Here Ca encodes the coupling to the slowly varying ID 
(scalar) potential field A{x), 



Ca^ — {d.^A/TTf + ipAA, 
m 

where pA denotes the total "gauge charge" density, 

PA = '^dxOaa/T^ ■ 

aa 

It is useful to define "charge" and "spin" boson fields, 
and "even" and "odd" flavor combinations. 



1 



± I 



'2m J 



(9) 



(10) 



(11) 



(12) 



with p. ~ p,a. Similar definitions hold for the (p fields. The 
commutation relations for the new 9, ip fields are unchanged. 
Integration over the gauge potential generates a mass term. 



CA=m{ep+~9'-°if 



(13) 



for the field 9p+ = J2aa ^qq/^- In the gauge theory analysis, 
we cannot determine the mean value 6*^*^, which is important 
for detailed properties of the SBM in Appendix A as well as 
for the discussion of nearby phases in Sees. IV B-IV E. But 
if we start with an interacting electron model, one can readily 
argue that the correct value in the SBM phase satisfies 
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(0) 



mod 27r 



(14) 



B. SBM by Bosonizing interacting electrons 



Consider then a model of electrons hopping on the zigzag 
strip. We assume that the electron hopping Hamiltonian is 
identical to the spinon mean field Hamiltonian, with first and 
second neighbor hopping strengths, ti,t2', 

H = - V[ii4(a:)c„(a; + 1) +<2ci(x)c„(a; + 2) +H.C.] 



X 



(15) 
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The electrons are taken to be at half-filling. The interactions 
between the electrons could be taken as a Hubbard repulsion, 
perhaps augmented with further neighbor interactions, but we 
do not need to specify the precise form for what follows. 

For t2 < 0.5ti, the electron Fermi sea has only one seg- 
ment spanning [— 7r/2, 7r/2], and at low energy the model is 
essentially the same as the ID Hubbard model. We know that 
in this case even an arbitrary weak repulsive interaction will 
induce an allowed four-fermion Umklapp term that will be 
marginally relevant driving the system into a ID Mott insula- 
tor. The residual spin sector will be described in terms of the 
Heisenberg chain, and is expected to be in the gapless Bethe- 
chain phase. 

On the other hand, for t2 > 0.5ti, the electron band has 
two Fermi seas as shown in Fig. 3. This is the case of pri- 
mary interest to us. As in the one-band case, Umklapp terms 
are required to drive the system into a Mott insulator But in 
this two-band case there are no allowed four-fermion Umk- 
lapp terms. While it is possible to study perturbatively the ef- 
fects of the momentum conserving four-fermion interactions 
and address whether or not the two-band metal is stable for 
some particular form of the lattice Coulomb repulsion, we do 
not pursue this here. Rather, we focus on the allowed eight- 
Fermion Umklapp term which takes the form. 

Ha = MCm]CRll^m]^R2l(^L^CLllCL2-\CL2l + H.c.) , 

(16) 

where we have introduced slowly varying electron fields for 
the two bands, at the right and left Fermi points. For repulsive 
electron interactions we have vs > 0. This Umklapp term is 
strongly irrelevant at weak coupling since its scaling dimen- 
sion is As = 4 (each electron field has scaling dimension 
1/2), much larger than the space-time dimension D ^ 2. 
To make progress we can Bosonize the electrons, just as we 



the latter case the electron excitations above the gap will 
correspond to instantons connecting adjacent minima of the 
cosine potential Eq. (17). In the spinon gauge theory there 
are no such fermionic excitations and the spinon excita- 
tions are linearly confined. This is appropriate for the spin 
model which has no "charge sector", and no notion of spinons. 
In the weak Mott insulating phase of the electron model, the 
Fermi wave vectors kpi, kp2 denote the momenta of the min- 
imum energy gapped electron excitations. What is the mean- 
ing, then, of the spinon Fermi wavevectors if the spinon ex- 
citations are unphysical? Within the spinon gauge theory the 
only gauge invariant (i.e. physical) momenta are the sums and 
differences of kpi, kp2, which correspond to momenta of the 
(low energy) spin excitations. In the electron model, the spin 
excitations below the charge localization length of the Mott 
insulator will be similar to that of electrons in the metal. On 
longer scales, the spin sector remains gapless, and this is the 
regime described below by the low energy effective theory 
of the SBM Mott insulator. It is these physical longer length 
scale spin excitations which are correctly captured by both the 
spinon gauge theory and interacting electron approaches. 



C. Fixed point theory of the SBM phase 



The low energy spin physics in either formulation can be 
obtained by integrating out the massive 9p+ field, as we now 
demonstrate. Performing this Gaussian integration leads to 
the effective fixed-point (quadratic) Lagrangian for the SBM 
spin liquid: 



^0 



SBM 



'-0 ' '-Q ' 



did for the spinons, cpaa 
Umklapp term becomes. 



+peaa) The eight-Fermion with the "charge" sector contribution. 



iJg = 2w8Cos(46lp+) 



(17) 



rP 



27rgo 



where as before 0p+ J^aa (^aaf^ and pe{x) = 2d^;9p+/Tr 
is now the physical slowly varying electron density. The 
Bosonized form of the non-interacting electron Hamiltonian 
is precisely the first part of Eq. (8), and one can readily con- 
firm that Ag = 4. But now imagine adding a strong density- 
density repulsion between the electrons. The slowly varying 
contributions, on scales somewhat larger than the lattice spac- 
ing, will take the simple form. Hp ^ Vpp1{x) ^ Vp{dxO p+Y . 
These forward scattering interactions will "stiffen" the Op^ 
field and will reduce the scaling dimension Ag. If Ag drops 
below 2 then the Umklapp term becomes relevant and will 
grow at long scales. This destabilizes the two-band metallic 
state, driving a Mott metal-insulator transition. The Opj^ field 
gets pinned in the minima of the potential, which gives 
Eq. (14). Expanding to quadratic order about the minimum 
gives a mass term of the form Eq. (13). For the low energy 
spin physics of primary interest this shows the equivalence 
between the direct Bosonization of the electron model and the 
spinon gauge theory approach. 

The difference between the spinon gauge theory and the in- 
teracting electron theory are manifest in the charge sector. In 



and the spin sector contribution. 



— {dreaaf+Va{dJa.f 



(18) 



(19) 



(20) 



The velocity vq in the "charge" sector depends on the product 
of the flavor velocities, vq = ^/viV2, while the dimensionless 
"conductance" depends on their ratio; 



90 



\/vi/v2 + \/v2/vi 



(21) 



Notice that go < 1, with go upon approaching the limit 
of a single Fermi surface {vi ^ Q,V2 ^ 0), and go 1 in the 
limit of two equally-sized Fermi surfaces {v2 /wi 1) that 
occurs when the two legs of the triangular strip decouple. 

In Sec. IV A, we also consider all symmetry allowed resid- 
ual short-range interactions between the low energy degrees 
of freedom and conclude that the above fixed-point theory can 
indeed describe a stable phase, with the only modification that 
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±2fcFi; 


±7r/2 


±(fcF2 — kpi); 
±(3A:fi + kF2) 






s 






2 ^ 4s 


1,1,9 
2 ^ 4s ^ 4 


1 


subd. 


B 


1 


1 


.9 


X 




subd. 








subd. 



TABLE I: Spin Bose-Metal fixed-point theory: Scaling dimensions 
of the spin S, bond energy B, and chirality x observables at various 
wavevectors Q in the top row. Entries with subdominant power laws 
are listed as "subd.". 



5o ~* 3 is now a general Luttinger parameter. Stability re- 
quires g < 1. There are also three marginal interactions that 
need to have appropriate signs to be marginally irrelevant. 

The gapless excitations in the SBM lead to power law cor- 
relations in various physical quantities at wavevectors con- 
necting the Fermi points. Here and in the numerical study 
Sec. Ill, we focus on the following observables: spin S{x), 
bond energy B{x) (i.e., VBS order parameter), and spin chi- 
rality xix): 

B{x) = S{x) ■ S{x + 1) , (22) 
Xix) = Six - 1) • [Six) X Six + 1)] . (23) 

In Appendix A, we give detailed expressions in the continuum 
theory. The most straightforward contributions are obtained 
by writing out, e.g.. Six) /t(a;)(T/(a;) in terms of the con- 
tinuum fermion fields and then bosonizing [see also Eqs. (A7) 
and (A9) for S(a;) andx(a;)]- We expect dominant power laws 
at wavevectors ±2kpa and ±7r/2 = Ti^pi + kp2), origi- 
nating from fermion bilinears composed of a particle and a 
hole moving in opposite directions. Such bilinears become 
enhanced upon projecting down into the spin sector (i.e. upon 
integrating out the massive 9p+ in the Bosonized field theory), 
and it is possible to compute the scaling dimension of any op- 
erator in terms of the single Luttinger parameter, g. It is also 
important to consider more general contributions, e.g., con- 
taining four fermion fields; this is best done using symmetry 
arguments and the corresponding expressions can be found in 
Appendix A. 

Table I summarizes such analysis of the observables by list- 
ing scaling dimensions at various wavevectors. We describe 
power law correlation of a given operator A at a wavevector 
Q by specifying the scaling dimension A^^ defined from the 
real-space decay 



(A(x)A(O)) ~ 



(24) 



The corresponding static structure factor (i.e. Fourier trans- 
form) has momentum-space singularity - ' ~ ^i2Aa„ -i 

The Q 
tions 



entries in Table I come from simple identifica- 



5'q=0 ^ dxdlcr + dx02a 

XQ=o ^ dx(fp~ ■ 



(25) 
(26) 
(27) 



In particular, the last line provides physical meaning to the 
"p— " sector - this spin-singlet sector encodes low energy fluc- 
tuations of the chirality. A direct way to observe the propagat- 
ing p— mode would be to measure the spectral function of the 
chirality, while in the present DMRG study we detect it by 
a |g| (i.e., V-shaped) behavior in the static structure factor at 
small wavevector q. 



III. DMRG STUDY OF THE SPIN BOSE-METAL IN THE 
Ji ./2 Aring MODEL ON THE ZIGZAG CHAIN 

We study the ring model Eq. (1) on the two-leg triangular 
strip shown in Fig. I. We use the ID chain picture and take 
site labels x ~ 1, . . . ,L where L is the length of the system. 
We use exact diagonalization (ED) and density matrix renor- 
malization group (DMRG)^*^"*^ methods supplemented with 
variational Monte Carlo (VMC)'''''*' to determine the nature 
of the ground state of the Hamiltonian Eq. (1). 



A. Measurement details 

We first describe numerical measurements. All calcula- 
tions use periodic boundary conditions in the x direction. In 
the ED, we can characterize states by a momentum quantum 
number k. On the other hand, our DMRG calculations are 
done with real-valued wavefunctions. This gives no ambigu- 
ity when the ground state carries momentum or tt and is 
unique. However, if the ground state carries nontrivial mo- 
mentum k ^ 0, TT, then its time-reversed partner carries —k, 
and the DMRG state is some combination of these. While the 
real space measurements depend on the specific combination 
in the finite system, the momentum space measurements de- 
scribed below do not depend on this and are unique. Most of 
the calculations are done in the sector with S^^^. = 0, which 
contains any ground state of the SU(2)-invariant system. 

The DMRG calculations keep more than m = 3200 states 
per block''^'"* to ensure accurate results, and the density ma- 
trix truncation error is of the order of 10^^. Typical relative 
error for the ground-state energy is of the order of 10"'' or 
smaller for the systems we have studied. Using ED, we have 
confirmed that all DMRG results are numerically exact when 
the system size is L = 24. The DMRG convergence depends 
strongly on the phase being studied, the system size, the type 
of the correlations, and the distance between operators. In 
the Bethe-chain and VBS-2 phases there is still good conver- 
gence for size L = 192, while in the SBM we are limited 
to i = 96 — 144 systems. The entanglement entropy calcu- 
lations are done with up to m = 6000 states in each block, 
which is necessary for capturing the long range entanglement 
in the SBM states where we find an effective "central charge" 
c~ 3. 

We have already specified the main observables in Sec. II C 
[cf. Eqs. (22-23)]. We measure spin correlations, bond en- 
ergy (dimer) correlations, and chirality correlations defined as 
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follows: 

C{x,x') = {S{x)-S{x')) , (28) 

D{x,x') = {B{x)B{x')) - {B)\ (29) 

X{x,x') = (x(a;)x(x')) ■ (30) 

For simplicity, we set D{x, x') = if bonds B{x) and B{x') 
have common sites, and similarly for X{x^ x'). Structure fac- 
tors C{q), D{q), wAX{q) are obtained through Fourier trans- 
formation: 



0(g) = -^^0(x,x')e- 



■iq{x—x ) 



(31) 



need the stability analysis of the SBM in Sec. IV, while the 
DMRG results are discussed afterwards in Sec. V. 

We explore more finely a cut i^^ring = Ji through the phase 
diagram Fig. 2, and our presentation points are taken from 
this cut. The Bethe-chain and VBS-2 phases are fairly con- 
ventional (for this idling, the VBS-2 is close to the decoupled- 
legs state at large J2 values). Nevertheless, it is useful to see 
measurements in these phases for comparisons. Such exam- 
ples are given in Appendix C, while here we focus on the Spin 
Bose-Metal point deep in the phase. We will discuss more 
difficult parts of the phase diagram Fig. 2 once we have the 
overall picture of the SBM. 



where O = C,D,X. We have 0{q) = 0{—q) and usually 
show only < g < tt. The spin structure factor at g = gives 
the total spin in the ground state: 



C{q = 0) 



{Si,) 



*S'tot(*S'tot + 1) 

L 



(32) 



In all figures, we loosely use "{BqB-qY' and "(xgX-g)" to 
denote D{q) and X{q) respectively. 

Turning to the VMC calculations, the trial wavefunctions 
are described in broad terms in Sec. I and in more detail in 
Appendix B. The states are labeled by occupation numbers of 
the two Fermi seas, (iVi, 7V2). Since iVi + iV2 = L/2, there 
is only one variational parameter There are three distinct 
regimes: i) Ni = L/2, N2 — 0, i.e., a single Fermi sea, which 
is appropriate for the Bethe-chain phase; ii) Ni ^ N2 ^ 
appropriate for the SBM; and Hi) Ni ~ N2 ~ L/A, i.e., de- 
coupled legs, which is a reasonable starting point for the large 
J 2 limit. 

In Appendix B, we describe correlations in the generic 
SBM wavefunctions and identify characteristic wavevectors 
2fc_Fa ~ 27rA^a/i, a = 1,2, and also 7r/2 (see Sec. II and 
Table I). One observation is that such wavefunctions corre- 
spond to a special case g = 1 in the SBM theory and thus can- 
not capture general exponents. Despite this shortcoming, the 
wavefunctions capture the locations of the singular wavevec- 
tors observed in the DMRG. We also try to improve the wave- 
functions by using a "gapless superconductor" modification 
described in Appendix B 2 and designed to preserve the sin- 
gular wavevectors while allowing more variational parame- 
ters. This indeed improves the trial energy and provides better 
match with the DMRG correlations at short scales, even if the 
long distance properties are still not captured fully. When pre- 
senting the DMRG structure factors, we also show the corre- 
sponding VMC results for wavefunctions determined by min- 
imizing the trial energy over the described family of states. 

Using the DMRG, we find four distinct quantum phases 
in the J2/J1 — A'^ing/Ji plane as illustrated in Fig. 2. In 
the small idling region, we have the conventional Bethe-chain 
phase at small J2 and Valence Bond Solid state with period 2 
(VBS-2) at larger J2- The SBM phase emerges in the regime 
idling > 0.2Ji and dominates the intermediate parameter 
space. Inside this region, we discover a new VBS state with 
period 3 (VBS-3). To fully understand the VBS-3 (in particu- 
lar, its relationship to the flanking Spin Bose-Metals) we will 



B. Representative Spin Bose-Metal points 

Proceeding along the A'ring = Ji = 1 cut through Fig. 2, 
we start in the Bethe-chain phase at large negative J2 (a repre- 
sentative point is discussed in Appendix C 1). As we change 
J2 towards positive value, the system undergoes a transition 
at J2 = —0.6. The new phase has characteristic spin correla- 
tions that are markedly different from the Bethe-chain phase. 
Figure 4 shows a representative point J2 = 0. The DMRG 
calculations are more difficult to converge and are done for 
smaller size L = 144 than in the Bethe-chain phase example 
(see also the entanglement entropy discussion below). 

Comparing with the Bethe-chain state (e.g.. Fig. 22 in Ap- 
pendix CI), there is no q = tt dominance in the spin structure 
factor. Instead, we see three singular wavevectors located at 
11 X 27r/144, 7r/2, and 61 x 27r/144. Our Gutzwiller wave- 
functions determined from the energetics have {Ni,N2) = 
(61, 11), and the corresponding 2kp2, —kp2 — kpi = 7r/2 
mod 27r, and 2kFi match precisely the DMRG singular 
wavevectors. The improved Gutzwiller wavefunction repro- 
duces crude short-distance features better, but it has the same 
long-distance properties as the bare Gutzwiller. As discussed 
in Appendix B, our Gutzwiller wavefunctions do not capture 
all power laws predicted in the general analytical theory. The 
wavefunctions appear to have equal exponents for spin corre- 
lations at these three wavevectors, while the theory summa- 
rized in Table I gives stronger singularities at 2fci?i, 2fci?2 and 
a weaker singularity at 7r/2. Very encouragingly, these the- 
oretical expectations are consistent with what we find in the 
DMRG, where we can visibly see the difference in the be- 
haviors at these wavevectors, particularly when we reference 
against the VMC results. 

Similar discussion applies to the bond energy (dimer) cor- 
relations shown in the middle panel of Fig. 4. The dominant 
features are at 2fci?i and 7r/2, and we also see a peak at a 
wavevector identified as AkF2, which is indeed expected from 
the SBM theory, cf. Table I. The theory predicts similar singu- 
larities at 2k PI and 2k p2, but for some reason we do not see 
the 2kp2 in the DMRG data, even though it is clearly present 
in the bare Gutzwiller. We suspect that this is caused by the 
narrowness in energy of the second Fermi sea when its popu- 
lation is low, so the amplitude of the 2fc^2 bond energy fea- 
ture may be much smaller The 2kp2 can still show up in 
the bare Gutzwiller since, as we discuss in Appendix B, this 
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FIG. 4: (Color online) Spin, dimer, and chirality structure factors at a 
representative point in the Spin Bose-Metal phase, A'ring = Ji = 1, 
J2 = (close to the Bethe-chain phase), measured in the DMRG 
for system size L — 144. We also show the structure factors in 
the Gutzwiller projection of two Fermi seas (A^i, N2) — (61, 11), 
and in the improved Gutzwiller wavefunction with parameters /o — 
1, fi — 0.65, /2 = —0.5 (see Appendix B; the parameters are de- 
termined by optimizing the trial energy within the given family of 
states). Vertical lines label important wavevectors expected in the 
SBM theory for such spinon Fermi sea volumes. We discuss the 
comparison of the DMRG, VMC, and analytical theory in the text. 



FIG. 5: (Color online) Spin, dimer, and chirality structure factors at a 
representative point in the Spin Bose-Metal phase, A'ring = Ji = 1, 
J2 = 3.2 (between the VBS-3 and VBS-2 phases) measured in 
the DMRG for system size L = 144. We also show results in 
the Gutzwiller projection of two Fermi seas (A^i, N2) — (44, 28), 
and in the improved Gutzwiller wavefunction with parameters fo = 
1, /i = 0.75, /2 = —0.1 (see Appendix B for details). Vertical lines 
label important wavevectors expected in the SBM theory. We dis- 
cuss the comparison of the DMRG, VMC, and analytical theory in 
the text. 
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wavefunction knows only about the Fermi sea sizes and not 
about the spinon energy scales like bandwidths or Fermi ve- 
locities. The improved Gutzwiller clearly tries to remedy this, 
although within its limitations. The Akp-z feature is not asso- 
ciated solely with the second Fermi sea and is less affected by 
this argument; indeed, 4fcF2 = — 4fcFi, and both Fermi seas 
"participate" in producing this feature as can be seen from 
Eq. (A31). 

Turning to the chirality structure factor in the bottom panel 
of Fig. 4, we see a feature at 7r/2, which in the SBM the- 
ory is expected to have the same singularity as the spin and 
dimer at this wavevector. We also see features at wavevec- 
tors kp2 — kpi and 'ikpi + kp2 in all our observables; these 
features are expected to be ~ \5q\ (i.e., V-shaped) in the 
Gutzwiller wavefunctions but have weaker singularity in the 
SBM theory, which is reasonably consistent with the DMRG 
measurements. Very notable in the chirality structure factor 
is a \q\ shape at small q. This can be viewed as a direct 
evidence for the gapless "p— " mode in the SBM, cf. Eq. (27). 
On the other hand, a feature at tt is weaker than V-shaped, in 
contrast with the Gutzwiller wavefunctions but in agreement 
with the SBM theory expectations in Table 1. 

To summarize, the correlations in the SBM phase are dra- 
matically different from the Bethe-chain phase, and we can 
match all the characteristic wavevectors using the Gutzwiller 
wavefunctions projecting two Fermi seas. We also understand 
the failure of the wavefunctions to reproduce the nature of 
the singularities and the amplitudes, while the Bosonization 
theory of the SBM is consistent with all DMRG observations 
even when the wavefunctions fail. 

With further increase of J2, continuing along the A'ling ~ 
Ji = 1 cut through Fig. 2, the SBM phase becomes promi- 
nently unstable towards the VBS-3 phase occupying the range 
1.5 < J2 < 2.5 and described in Sec. V A. Interestingly, as 
we increase J2 still further, the SBM phase reappears with its 
characteristic correlations shown in Fig. 5 for a representative 
point J2 = 3.2. Much of the SBM physics discussion that we 
have just done at J2 = carries over here, with the appropri- 
ately shifted locations of the singular wavevectors. The singu- 
larities at wavevectors qiow = 2fci?2 andqhigh = 2kpi are now 
more equally developed and are detected in the spin as well 
as dimer structure factors. The two wavevectors are closer 
to 7r/2 and are located symmetrically in accordance with the 
general "sum rule" 2kpi + 2kp2 = tt, while the compara- 
ble strengths reflect the more similar energy bandwidths of 
the two Fermi seas. The apparent lack of the wavevector tt/2 
in the DMRG dimer structure factor is similar to that in the 
Gutzwiller wavefunction and is a matrix element effect for the 
first-neighbor bond when the sizes of the two Fermi seas ap- 
proach each other On the other hand, the strength of the 4fci?2 
dimer feature is very notable here; it can indeed be dominant 
in the SBM theory for sufficiently small Luttinger parameter 
g, cf. Table I. The improved Gutzwiller wavefunction modi- 
fies the structure factors in the right direction but clearly does 
not succeed reproducing them accurately - as noted before, 
our wavefunctions do not contain the full physics expected in 
the Bosonized theory. 

One technical remark that we want to make here is that 



the DMRG ground state at this point J2 = 3.2 and size 
L = 144 appears to have non-trivial momentum quantum 
number fc ^ 0, tt. We deduce this by observing that the mea- 
sured correlations 0(.t, x') depend not just on x — x' but on 
both X and x' , and by seeing characteristic beatings as a func- 
tion of X and x' (while the (/-space structure factors are well- 
converged). On the other hand, the VMC wavefunction shown 
in Fig. 5 has momentum zero (see Appendix B) and all mea- 
surements depend on a; — x' only. If we assume that the beat- 
ings originate from the DMRG state being a superposition of 
|fc) and I — fc), we can extract 2fc and find it to be consistent 
with the state |fc) constructed from the VMC by moving one 
spinon across one of the Fermi seas (2fc ~ Akp2 = —Akpi). 
It is plausible that such state happens to have a slightly lower 
energy in the given finite system (e.g., at the same J2 = 3.2, 
we find trivial k for L ~ 72 but non-trivial k for L = 96, 
likely reflecting finite-size effects on the filling of the last few 
spinon orbitals). We have not attempted to construct a trial 
spin-singlet state with the right momentum quantum number 
for the present L. Still, we expect that the structure factors are 
not very sensitive to such rearrangements of few spinons in the 
large system limit. Indeed, we find that the structure factors 
have the same features for different system lengths L = 72, 
96, and 144. It is also worth repeating that our structure factor 
measurements using Eq. (31) do not depend on which specific 
combination of \k) and | — A:) is found by the DMRG proce- 
dure. 



C. Evolution of the singular wavevectors in tlie SBM 

We further illustrate the Spin Bose-Metal by showing evo- 
lution of the DMRG structure factors and singular wavevec- 
tors as we move inside the phase. The spin and dimer struc- 
ture factors are shown in Fig. 6 for L = 96 and varying pa- 
rameter J2/J1 inside the SBM phase adjacent to the Bethe- 
chain phase. With increasing J2, the singular wavevector 
9higii (identified as 2kpi in the VMC) is moving away from 
TT towards smaller values, while the singular wavevector giow 
(identified as 2kp2) is moving to larger values; this corre- 
sponds to spinons being transfered from the first to the second 
Fermi sea as found in the VMC energetics. The spin and dimer 
correlations show similar behavior at 2fci?i , and both have fea- 
tures also at the wavevector tt/2. The lack of visible dimer 
feature at 2kp2 was discussed for the point J2 = earlier, 
and it is likely that something similar is at play here. On one 
hand, the second band is narrow when we just enter the SBM 
since the second Fermi sea starts small. On the other hand, 
in the region labeled "partial FM"close to the VBS-3 phase in 
Fig. 2, the DMRG finds nonzero magnetization in the ground 
state. We think that this occurs in the second Fermi sea and 
indicates an effective narrowness of this band near the VBS-3 
phase as well, so the 2k p2 energy feature may indeed be weak 
in the whole SBM phase between the Bethe-chain and VBS-3. 

We think that the same physics also starts causing conver- 
gence difficulties in the DMRG for the L = 96 systems shown 
in Fig. 6. Specifically, we can use Eq. (32) to extract the 
ground state spin Stot and find non-integer values of order 
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FIG. 6: (Color online) Evolution of the structure factors in the Spin 
Bose-Metal phase between the Bethe-chain and VBS-3 (or partial 
FM), measured by the DMRG for system size L — 96. We can 
track singular wavevectors ("Bose surfaces") as spinons are trans- 
fered from the first to the second Fermi sea upon increasing J2 . In 
the spin structure factor, the wavevector (jhigh that starts near n is 
identified as 2fcFi and the wavevector giow that starts near is iden- 
tified as 2fcF2 (see Fig. 4 and text for details). The ghigh and giow are 
summarized in Fig. 8. 



1, e.g., 5'tot(5'tot + 1) = 0.3, 1.0, 1.9, 2.1, 3.4 for the points 
J2 0, 0.2, 0.4, 0.6, 0.8 in Fig. 6. This can only happen if the 
DMRG is not converged and is mixing several states that are 
close in energy but have different spin. Correspondingly, we 
observe a difference between {S^S tg) and {S^S^g) structure 
factors (recall that we are working in the sector Sf^j ~ 0). We 
do not show these graphs, but the difference is localized near 
2kp2, where the (SgStq) has much sharper feature while the 
{SgS^g) has weaker feature. On the other hand, there is es- 
sentially no difference near 2kFi- This strongly suggests that 
the origin of the convergence difficulties lies with the second 
Fermi sea. The {Sq ■ structure factor that is shown in 
Fig. 6 does not mix different ^tot states but only sums the 
corresponding structure factors and is less sensitive to these 
convergence issues. The fact that the peaks in the top panel of 
Fig. 6 are located symmetrically with respect to 7r/2 suggests 



FIG. 7: (Color online) Evolution of the structure factors in the Spin 
Bose-Metal phase between VBS-3 and VBS-2, measured by the 
DMRG for system size L — 96. The ghigh = ^kpi and giow ~ 2A:f2 
bracket the 7r/2 and approach each other with increasing J2 (they are 
summarized in Fig. 8). Very notable here is the strong 4A;f2 peak in 
the dimer structure factors that evolves out of the Bragg peak present 
in the VBS-3 phase at 27r/3 (see Sec. VA and Fig. 15); the 2A:fi 
moves away from the 27r/3 in the opposite direction. 

that this region is still spin-singlet SBM (but is on the verge of 
some magnetism in the second band). Finally, the ^tot values 
for the same parameters but system size L = 48 are indeed 
well converged to zero, however with increasing number of 
states needed in the DMRG blocks for larger J2. We thus 
conclude that the points in Fig. 6 are spin-singlet SBM. The 
DMRG convergence difficulties for the larger L are in accord 
with the presence of many low-energy excitations (see also 
the discussion of the entanglement entropy below). We will 
further test our intuition that this region is close to some weak 
subband ferromagnetism in Sees. V B and VI by adding anti- 
ferromagnetic or ferromagnetic J3 to suppress or enhance the 
FM tendencies. 

Consider now Fig. 7 that shows evolution of the structure 
factors in the SBM phase between the VBS-3 and VBS-2. 
The Qhigh = 2kpi and giow = 2kp2 continue moving to- 
wards each other with increasing J2, and the spin structure 
factors become nearly symmetric with respect to 7r/2. When 
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the ghigh and q\ow peaks merge at 7r/2, which in the VMC 
would correspond to decoupled legs, one expects''^ ''' that a 
new instability will likely emerge resulting in a VBS-2 state 
(we discuss a representative point in Appendix C 2). 

A very notable feature in the SBM dimer structure factor is 
the strong 4:kF2 peak. Foretelling a bit, this peak can be traced 
as evolving out of the dimer Bragg peak at 27r/3 of the VBS-3 
phase to be discussed in Sec. V A. Turning this around and 
approaching the VBS-3 phase by decreasing J2, we can view 
the VBS-3 as an SBM instability when the dimer Akp2 peak 
merges with the 2kpi singularity, 4fcj?2 = 2kpi = 27r/3. 

Finally, in this SBM region the DMRG converges well to 
spin-singlet states. The remark we made when discussing the 
point J2 = 3.2 applies to all points shown in Fig. 7 with L = 
96; they show correlations 0{x, x') beating in both x and x', 
which can be interpreted similarly to the earlier J2 = 3.2 case 
by assuming superposition of degenerate ground states with 
opposite momentum quantum numbers. 

Fig. 8 summarizes the singular spin wavevectors extracted 
from plots Uke Figs. 6 and 7, superimposed on the phase di- 
agram of the model along the cut /^ling = Ji- Remark- 
ably, the singular wavevectors throughout the entire SBM 
phase are well captured by the improved Gutzwiller wave- 
functions, as we have illustrated in Figs. 4 and 5. These singu- 
lar wavevectors are intimately connected to the sign structure 
of the ground state wavefunction, indicating a striking coinci- 
dence between the exact DMRG ground state wavefunction 
and the Gutzwiller projected VMC wavefunction. Besides 
the SBM regions. Fig. 8 also shows the Bethe-chain (cf. Ap- 
pendix C 1), VBS-2 (Appendix C 2), and VBS-3 (Sec. VA) 
phases. 

We now mention more difficult points in the overall phase 
diagram. The lightly hatched SBM region in Fig. 8 indicates 
the discussed rising DMRG difficulties of not converging to 
an exact singlet for L = 96. Such DMRG states are shown as 
open circles with crosses in Fig. 2. As we have already men- 
tioned, the estimated S'tot values are not converged and are of 
order 1 for L = 96 (but are converged to zero for L = 48), 
while (/low and Qhigh are located symmetrically around tt/2; 
all this suggests that the phase is spin-singlet SBM. On the 
other hand, at points J2 = 1.2—1.5 not marked in Fig. 8 
but shown as star symbols in Fig. 2, the estimated S'tot val- 
ues are larger and the apparent dominant wavevectors are no 
longer located symmetrically. Here we suspect a modification 
of the ground state, likely towards partial polarization of the 
second Fermi sea; this "partial FM" region is also indicated 
by cross-hatching in Fig. 8. 



D. Entanglement entropy and effective central charge of the 
SBM 



We explore properties of the SBM phase that can further 
distinguish it from the Bethe-chain and VBS states. Earlier 
we have noted that we need to keep more states per block to 
achieve similar convergence for the SBM phase in compari- 
son with the Bethe-chain and VBS phases, indicating stronger 
entanglement between subsystems in the SBM. Bosonization 
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FIG. 8: (Color online) Cut through the phase diagram Fig. 2 at 
A'ring/Ji = 1, showing evolution of the most prominent wavevec- 
tors in the spin structure factor. In the Bethe-chain phase we have 
singular antiferromagnetic qaf ~ ti", cf. Fig. 22. In the Spin Bose- 
Metal we have singular ghigh ~ ^kpi and giow ~ 2A:f2 located 
symmetrically about 7r/2, cf. Figs. 6 and 7. In the VBS-3 we have 
singular qaf ~ tt/S, cf. Fig. 15. In the VBS-2 region for large 
J2, cf. Fig. 23, we have dominant correlations at 7r/2 corresponding 
to the decoupled legs fixed point, which is likely to be unstable to- 
wards opening a spin gap."* " The dotted lines show results for the 
improved Gutzwiller wavefunction. 



analysis finds that the SBM fixed point theory has three free 
Boson modes. One can associate a central charge 1 with each 
mode. Despite the fact that they have different velocities (so 
the full system is not conformally invariant), we expect that 
the total entanglement entropy should have a universal behav- 
ior described by a combined central charge c — 3. 

In general, for a one-dimensional gapless state with confor- 
mally invariant correlation functions in space-time, the entan- 
glement entropy for a finite subsystem of length X inside a 
system of length L with periodic boundary conditions varies 
as« 



S{X,L) = ^ log (^sin^ 

3 \TT L 



(33) 



where ^ is a constant (independent of the subsystem length) 
and c is the effective central charge. The virtue of the entan- 
glement entropy is that it does not depend on the mode veloc- 
ities and in principle measures the number of gapless modes 
directly from the ground state wavefunction. 

Figure 9 shows the entanglement entropy S{X, L) as a 
function of X for different quantum phases for a finite system 
length L = 72. The results are obtained from the DMRG for 
representative points taken from the same cut Kiing ~ Ji ~ 1 
discussed earlier The entropy can be well fitted by the ansatz 
Eq. (33) with different c values. 

The Bethe-chain state (at J2 = —1) gives central charge 
c = 1.0 consistent with one gapless mode. On the other hand, 
the entanglement entropy for either of the two SBM examples 
J2 = and J2 = 3.2 is much larger and can be fitted by close 
values c = 3.1 and c = 3.2, respectively. The closeness of the 
central charges in these two different SBM states (cf. Figs. 4 
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FIG. 9: (Color online) Entanglement entropy at representative points 
in the Bethe-chain (J2 = -1), SBM (J2 = and J2 = 3.2), 
VBS-3 (J2 = 2), and VBS-2 (J2 = 4) phases, talcen from the 
cut A'ring — Ji = 1, measured in the DMRG for system size 
L = 72 with periodic boundary conditions. We use Eq. (33) to fit 
data over the range 6 < X < 66, which gives the best estimates as 
c = 1.0 (Bethe-chain), c = 3.1 (SBM at J2 = 0), c = 3.2 (SBM at 
J2 = 3.2), c = 1.6 (VBS-3), and c = 2.1 (VBS-2). 
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FIG. 10: (Color online) Entanglement entropy for various system 
sizes for the Bethe-chain (J2 = —1) and SBM (J2 = 0) points, 
cf. Fig. 9, plotted versus scaling variable rf = ^ sin We also 
show fits to Eq. (33) done for the larger sizes. The Bethe-chain data 
collapses well and is fitted with c = \. The SBM data is fitted with 
c = 3.1 (we show the same fit as in Fig. 9 for this SBM point). The 
collapse of different L is less good, which is likely due to impre- 
cise scaling of the discrete shell filling numbers with L (see text for 
details). 



and 5) indicates the universal behavior of the entanglement 
which is independent of the details like the relative sizes of 
the spinon Fermi seas. 

Interestingly, the VBS-2 point at J2 = 4 is fitted by c = 2.1, 
which is related to the fact that the wavefunction is close to the 
decoupled-legs limit (see Appendix C 2 and Fig. 23), and this 
L = 12 system "does not know" about the eventual spin gap 
and very small dimerization. 

Finally, the fitted effective central charge for the VBS-3 ex- 
ample is around c = 1.6. The oscillatory behavior of SiX) 
reflects translational symmetry breaking in the DMRG state. 
Note that the entropy values are larger here compared with 
the Bethe-chain or VBS-2 cases, which is probably due to a 
mix of degenerate states in the DMRG wavefunction. How- 
ever, the overall X-dependence is clearly weaker than in the 
VBS-2 and is approaching the Bethe-chain behavior for large 
X. 

To better understand finite size effects, we focus on the 
SBM and Bethe-chain phases and discuss the universal de- 
pendence of the entropy on the scahng variable d = ^ sin . 
Figure 10 shows L) as a function of d for several sys- 
tem sizes for the SBM point J2 = and the Bethe-chain point 
J2 = -1. 

At the SBM point, the data for the two larger sizes L = 60 
and 72 collapse onto one curve, which can be reasonably fit- 
ted by S{X,L) = (3.1/3) log((i) + 0.88, strongly suggest- 
ing the effective central charge c ~ 3. The smaller sizes 



L = 24 and 48 have somewhat shifted entropy values com- 
pared with the L = 60 and 72 collapse, but show roughly 
similar slope for the largest d. The differences are likely due 
to finite-size shell filling effects. Indeed, we can measure the 
structure factors and characterize the presumed SBM states by 
the spinon occupation numbers of the two Fermi seas; we find 
{Ni,N2) = (10,2), (20,4), (26,4), and (31,5) fori = 24, 
48, 60, and 72 respectively, and these numbers in each Fermi 
sea do not precisely "scale" with L. 

On the other hand, the Bethe-chain case does not have such 
effects and data for all sizes collapse. The results can be well 
fitted by 5(X,i) = (1/3) log((i) + 0.99 as shown in the same 
figure for system sizes i = 24 to 96. We note that while 
the entropy for the Bethe-chain phase for L = 96 is fully 
converged by keeping up to m = 4200 states in the DMRG 
block, the entropy for the SBM for L = 72 is still increasing 
slowly with the number of states kept, and we estimate that the 
error in S{X = L/2,L) is around a few percent when m = 
6000 states are kept per block (comparing to an extrapolation 
to m = 00). Indeed, the SBM data for L = 72 is bending 
down slightly from the fitted line at the larger d corresponding 
to X ~ L/2, as can be seen in Fig 10, which is probably 
because the data is less converged. 

To summarize, the entanglement entropy calculations es- 
tablish the SBM as a critical phase with three gapless modes 
and clearly distinguish it from the Bethe-chain and VBS-2 



14 



phases (or the decoupled-legs limit). We also note that the 
structure factor measurements and detection of all features as 
in Figs. 6 and 7 did not require as much effort and was done 
for larger sizes than the entropy; however, to characterize the 
long distance power laws accurately one would probably need 
to capture all entanglement, which we have not attempted. 



IV. STABILITY OF THE SPIN BOSE-METAL PHASE ; 
NEARBY PHASES OUT OF THE SBM 

A. Residual interactions and stability of the SBM 

We account for residual interactions between low-energy 
degrees of freedom in the SBM theory. Sec. 11, by considering 
all allowed short-range interactions of the spinons. The four- 
fermion interactions can be conveniently expressed in terms 
of chiral currents: 



The remaining exponentials only depend on the fields Baa 
and ipacr so that the charge and spin sectors decouple. Since 
the rest of is quadratic, C ~ Cq + takes the same form 
as Cq in Eq. (19) except with go,vo — > g,v, where 

2 ^ ivi-KiM{v2-K2M-{K,M' ,,,, 

^ {v^+V2^2Xl,/7Ty~{XlJn + X^,,/ny'' ' 

V = ^-{vi+V2 + \1jTr + KjTi-2\1jTi) . (42) 

In the spin sector, the remaining interactions are given by, 

-^1 = — AJ^„ JfJaa ■ J Laa — Xi2{J Rll ' «/l22 + Jill • J R22) ■ 
a 

(43) 

When we write this in the bosonization, the J£ pieces con- 
tribute to the harmonic part of the action 



A 



Jpab — fpaafPba '7 Jpab ~ -^.fpaa'^afS fpbp 



(44) 



(34) 



We assume that interactions that are chiral, say involving only 
right movers, can be neglected apart from velocity renormal- 
izations. The most general four-fermion interactions which 
mix right and left movers can be succinctly written as, 

= ^ (w^l^jRabJLab + ^'abJ RaaJ Lbb) , (35) 
a,b 



+ ^mH^la){d^V2a) ~ {dJla){dj2a)] ,(45) 

while the J^J]^ + J^ J^ produce nonlinear potential 

V± = Y.^aaCOs{2V20aa) (46) 

a 

+ 2A?2rcos(26l^+)cos(2(^^_) . (47) 



£J ^ - ^ [w:Jnab ■ JLab + KbJRaa ' JLbb) ,(36) ^ 1-loop RG analysis gives the following flow equations. 



with wii = W22 = (convention), wi2 = W21 (from Her- 
miticity), and A12 = A21 (from L symmetry). There are 

8 independent couplings: ^^'2'^' ^11'^' '^22'^' ^^\2 ■ 

We treat these interactions perturbatively as follows. First 
we bosonize the interactions and obtain terms quadratic in 
dxdaa ™d dx^Paa ^s Well as tcrms involving products of four 
exponentials e^'^"° and e*"^"" . We next impose the condi- 
tion that 6p+ is pinned and compute the scaling dimensions of 
the exponential operators. 



iKaf 



dXl2 

di 



7r(wi + V2) 



(48) 



The w'^2" terms give 



W = {WI2JR12JL12 - W12JR12 ■ Jlu) + H.c. = 
= cos(2</jp_)< iwi2 cos(2(/3ct-) — rcos(20cr-) 



(37) 



(38) 



w^2 cos(2v?^_ ) + r cos(26l^_ ) + 2r cos(20^h 



.(39) 



where 



(40) 



The Wi2'^ terms have scaling dimension 1 + g^^ > 2 and are 
irrelevant in the bare theory, Eq. (18), and henceforth dropped. 
The detailed W expression will be used later when we analyze 
phases neighboring the SBM. 



When Xii, X22, and A'^2 are positive, they scale to zero and 
the quadratic SBM Lagrangian Cq^^, Eq. (18), is stable. We 
also require that the renormalized g is smaller than 1 so that 
the ^12'^ terms in Eq. (37) remain irrelevant. In Sec. IV B, 
we consider what happens when g > 1 or when some of the 
Aii,A22, and X12 change sign and become marginally rele- 
vant. 

The above stability considerations are complete for generic 
incommensurate Fermi wavevectors. At special commensu- 
rations, new interactions may be allowed and can potentially 
destabilize the SBM. Such situations need to be analyzed sep- 
arately, and in Sees. IV D and IV E we consider cases relevant 
for the VBS-3 and ChiraUty-4 phases found by the DMRG in 
the ring model. Sees. V A and V B. 

We want to make one remark about the allowed interac- 
tions, which will be useful later. Let us ignore for a mo- 
ment the pinning of 9p+; for example, let us think about 
the electron interactions in the approach of Sec. II B. Three 
of the eight 9 and (p fields, namely fp-i-, fa+, and 6'p_, 
do not appear as arguments of the cosines, and the action 
has continuous symmetries under independent shifts of these. 
The first two symmetries correspond to microscopic conser- 
vation laws for the total charge Q ^ J dx0p+ and the to- 
tal spin S'^ ^ J dxda+- On the other hand, the invariance 
under the shifts of 9p^ corresponds to the conservation of 
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X = Nni^ Nli- Nr2 +Nl2 - J d^ifp-, where Npa 
denotes the total number of fermions near Fermi point Pkpa- 
This is not a microscopic symmetry, but emerges in the con- 
tinuum theory for generic kpi, kp2- Indeed, writing the total 
momentum P = {Nri - NLi)kFi + {Nr2 - ^L2)fc_F2 = 
Xkpi — {Nii2 — AfL2)7r/2, we see that any attempt to change 
X violates the momentum conservation except for special 
commensurate kpi. We thus conclude that Lpp+ and tp„^ 
can never be pinned by the interactions, while 6'p+ can not 
be pinned generically except at special commensurate points. 



B. Gapped paramagnets when ,g > 1 

We now consider phases that can emerge as some instabili- 
ties of the Spin Bose-Metal. We use heavily the Bosonization 
expressions for various observables given in Appendix A. As 
we have already mentioned, when .g > 1 the interactions in 
Eq. (37) are relevant and the SBM phase is unstable. We can 
safely expect that as a result of the runaway flows, the vari- 
ables ifp- and 9a+ will be pinned. The situation is less clear 
with the remaining parts of the potential since we cannot pin 
simultaneously O^r- and ipa-. Still, it is possible that the situ- 
ation is resolved by pinning one variable or the other For ex- 
ample, depending on whether ^i^^ have the same or 
opposite signs, it is advantageous to pin O^r- or i^a-- If either 
pinning scenario happens, there remain no gapless modes in 
the system. It is readily established in the cases below that all 
spin correlations are short-ranged, i.e., we have fully-gapped 
paramagnets; also, (S^) 7^ in all cases, so the translational 
symmetry is necessarily broken. 

1. TOi2'U'i2 > and pinned ipp-, Ocr+, 6cr--' period-2 VBS 

Consider the case when the 9cr- is pinned. Using Ap- 
pendix A, we see that B-^ obtains an expectation value, while 
Bt^/2 and Xn/2 ^re short-ranged. It is natural to identify this 
phase as a period-2 Valence Bond Solid shown in Fig. 1 1 . The 
pinning values and therefore some details of the state will dif- 
fer depending on the sign of the coupling w^2' but in either 
case the ground state is two-fold degenerate. [Here and be- 
low, when we find pinning values of appropriate cp-s and 9-s 
minimizing a given potential, we determine which solutions 
are physically distinct by checking if they produce distinct 
phases modulo 27r in the bosonization Eq. (5). More prac- 
tically, following Ref. 31 Sec. IV.E.l, the chiral fermion fields 
remain unchanged under (paa ^fiaa + T^{(-Raa + f-Laa), 
daa daa + '^{(■Raa " (-Lao), whcrc ipaa Can be arbitrary 
integers. This gives redundancy transformations for the pi, 
CT± fields that we use to check if the minimizing solutions are 
physically distinct.] 

2. W12WI2 < and pinned ipp- , ^(t+, ^ct-.' period-4 structures 

Consider now the case when the (/s^- is pinned. We find 
that either e7r/2 in Eq. (A42) or x-k/2 in Eq- (A43), but not 




FIG. 11: Valence Bond Solid with period 2, where thicker lines in- 
dicate stronger bonds. To emphasize the symmetries of the state, we 
also show second- and third-neighbor bond energies, but details can 
be different in different regimes. For example, the VBS-2 state in 
the A'ring = case has dominant first-neighbor dimerization. On 
the other hand, in the model with A'ring = Ji = 1, Jz = 0.5, the 
putative VBS-2 region between the Bethe-chain and SBM phases in 
Fig. 16 has significant third-neighbor modulation but only very small 
first-neighbor one. 

both, obtains an expectation value. Thus we either have a 
period-4 VBS or a period-4 structure in the chiralities. Which 
one is realized depends on details of the pinning. 

As described in Appendix A, we work with the +1 eigen- 
state of the operator F [our Eqs. (A42-A43) already assume 
this]. With this choice, to minimize the potential in Eq. (37) 
we require 

cos(2v7^_) cos(26l^+) = ±1 . (49) 
Depending on the sign of w^j' we have: 

a) w;^2 > : cos(2v3p_) = - cos(26'^+) , (50) 

b) < : cos(2v3p_) = 008(26*^+ ) . (51) 

a) In this case, {ejr/2) ^ 0, {x-k/2) ~ 0, i.e., we find period- 
4 valence bond order Note that (e7r/2) can take four indepen- 
dent values (e7r/2) = e*" = e**'^/'*, ^^^'^^/^^ where we have 
assumed that Op+ is fixed by Eq. (14). To visualize the state, 
we examine the corresponding contributions to the first- and 
second-neighbor bond energies: 

5B^^x+i - cos(|x+ J + a) =={+,0,-,0,...}, (52) 

SBx,x+2 COS - + aj ={+,-,-,+,...}. (53) 

One can either use symmetry arguments or write out the mi- 
croscopic hopping energies explicitly to fix the phases as 
above [see Eq. (A8), which generalizes to n-th neighbor bond 
as Sq'^ ~ e^'^/^eg for Q ^ n]. Each line also shows 
schematically the sequence of bonds starting at a; = for 
a = — 7r/4. The four independent values of a correspond to 
four translations of the same VBS state along x. The pattern 
of bonds is shown in Fig. 12, where the more negative energy 
is associated with the stronger dimerization. When viewed on 
the two-leg ladder, this state can be connected to a state with 
independent spontaneous dimerization in each leg. 

b) Here we have (e^/2) = 0, {Xit/2) 0' i-^-, period-4 
structure in the chirality x- The pattern is 

X(x) - cos (^^x + a) ={+,-,-,+,...}. (54) 

There are four independent values of (Xt/2) = e*" = 
g±t7r/4^ g±j37r/4^ Corresponding to four possible ways to regis- 
ter this pattern on the chain. The state is illustrated in Fig. 13. 
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FIG. 12: Top: Valence Bond Solid with period 4 suggested as one 
of the instabilities out of the Spin Bose-Metal phase, Sec. IVB2a. 
Thick lines indicate stronger bonds. Bottom: In the two-leg triangu- 
lar ladder drawing, we see roughly independent spontaneous dimer- 
ization in each leg. 
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FIG. 13: Top: Chirality order with period 4 suggested as one of the 
instabilities out of the SBM phase, Sec. IV B 2b. In the ID chain pic- 
ture, the chirality pattern is given by Eq. (54); x(a^) is associated with 
the [x — 1 , a: , a; + 1] loop and the arrows on the links show one "gauge 
choice" to produce such "fluxes" in the spinon hopping. Bottom: In 
the two-leg triangular ladder drawing, we see alternating chiralities 
on the up-triangles along the strip and alternating chiralities on the 
down-triangles. 



When drawn on the two-leg ladder, chiralities on the upwards 
pointing triangles alternate along the strip, and so do chirali- 
ties on the downwards pointing triangles. 



3. Gapped phases in the spinon language 

With an eye towards what might happen in the 2D spin liq- 
uid, it is instructive to discuss the above phases in terms of the 
spinons. To this end, we can rewrite the W term, Eq. (37), as 
follows: 

W = K2+<2/4)[P;P2+H.c.] (55) 

- i<2 - y^l2/m.fLfL2)ifljB2) + H.C.] . (56) 

Here = /Lt/Lx " fLjL]' creates a "Cooper pair" 
in band a. The preceding two sections can be then viewed 
as follows. When u'^2'"^i2 > 0' we minimize the first line 
by "pairing and condensing" the spinons; once everything is 
done, we get the period-2 VBS state. On the other hand, when 
w 12^12 < 0, we minimize the second line by developing ex- 
pectation values in the particle-hole channel. Using 

(/ii/L2)(/li/fl2) + H.c. = 2[el^2^^/2-xl/2X./2] , (57) 

we get either the period-4 dimer or period-4 chirality order 
depending on the sign of ?i;^2- 



C. Nearby phases obtained when some of the • Jl 
interactions, Eq. (43), become marginally relevant 

Let us now assume g < 1, so the singlet "p— " sector is not 
a priori gapped. We consider what happens when some of the 
couplings Xii , X22 , A12 in Eq. (43) change sign and become 
marginally relevant. We analyze this as follows. Consider the 
potential V±, Eq. (47), again working with the +1 eigenstate 
of the operator F. If one (or several) of the couplings becomes 
negative, we have runaway flows Eq. (48) to still more nega- 
tive values. We then consider pinned field configurations that 
minimize the relevant part of the V± — this is what happens 
in the spin sector. Next we need to include the interactions 
Eq. (37) with the singlet "p— " sector, since they can become 
relevant once some of the "cr" fields are pinned. We now con- 
sider different possibilities. 



1. AJi > 0, A^2 > 0, A?2 < 

In this case, only the XI2 is relevant and flows to large neg- 
ative values. We therefore pin the fields 9^+ and (pa-- To 
minimize V±, the pinned values need to satisfy Eq. (49). The 
spin sector is gapped and all spin correlations decay exponen- 
tially; we also have (S^r) ^ 0, so the translational symmetry 
is broken. 

Next we include the interactions Eq. (37). Using Eq. (49), 
the important part is 

W = (4w^2 - ^^12) cos(26l^+) cos(2v3p_) . (58) 

The ipp^ is dynamical at this stage, but the 9c+ is pinned and 
the W now has scaling dimension 1/g. The possibilities are: 

a) g < 1/2: The W term is irrelevant and the singlet sector 
remains gapless. One manifestation of the gaplessness is that 
B±Tr/2 and x±7r/2 havc power law correlations characterized 
by scaling dimension l/(4g). Thus we have a coexistence of 
the static period-2 VBS order and power law VBS and chiral- 
ity correlations at the wavevectors ±7r/2. 

h) g > 1/2: The W term is relevant and pins the field ipp^ 
leaving no gapless modes in the system. Such fully-gapped 
situation has already been discussed in Sec. IV B 2. This gives 
either the period-4 VBS or period-4 chirality phase. 



2. AJi < 0, A^2 < 0, X12 > 

In this case, the X^i and A22 are relevant and flow to large 
negative values while A j2 is irrelevant. Then both di^ and 92a 
are pinned and satisfy cos(2^/2^^lCT) = cos(2V26'2cr) = 1. 
The spin sector is gapped and all spin correlations are short- 
ranged. All correlations at tt/2 are also short-ranged. The 
translational symmetry is broken since {6^^) 7^ 0. Including 
the interactions with the singlet sector as in the previous sec- 
tion, we have: 

a) If g < 1/2, the "p— " sector remains gapless and B2kFi 
and B2kF2 have power law correlations with scaling dimen- 
sion g/4. These coexist with the static period-2 VBS order 
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b) If 5 > 1/2, we also pin Lpp- and the situation is essen- 
tially the same as in Sec. IV B 1 . This gives the fully-gapped 
period-2 VBS phase. 



3. A?i < 0, \22 > 0, A?2 > 

In this case, only the AJj is relevant and pins die,. Spin cor- 
relations at 2fcFi and all correlations at 7r/2 are short-ranged. 

We now include the interactions Eq. (37); the importantpart 

is 

W = -{■iw^2 + 3^12) cos(V26'i^) 008(^26*2^) cos(2(/3p_) . 

(59) 

Both the "2fT" and modes are dynamical at this stage, 
and the W has scaling dimension 1/2 + I/5. 

a) g < 2/3: The W term is irrelevant and we have two 
gapless modes in this phase. (S, B)2kF2 have the same scaling 
dimension 1/2 + .g/4 as in the SBM phase, while B2kFi has 
scaling dimension g/4. Furthermore, (S", S)^ have scaling 
dimension 1/2. 

h) g > 2/3: The W term is relevant pinning both 92a and 
(fip- . This is the already encountered fully gapped period-2 
VBS state. 

The case with Xfi > 0, A22 < 0, AJ2 > is considered 
similarly. 

Finally, in the case \i2 < and either A^i < or A22 < 0, 
we can not easily minimize the potential Eq. (47) since we 
have non-commuting variables under the relevant cosines. We 
do not know what happens here, although one guess would be 
that one of the relevant terms wins over the others and the 
situation is reduced to the already considered cases. 

To summarize, we have found several phases that can be ob- 
tained out of the Spin Bose-Metal: 1) fully gapped period-2 
VBS; 2), 3) fully gapped period-4 phases, one with bond en- 
ergy pattern and the other with chirality pattern; 4) period-2 
VBS coexisting with one gapless mode in the singlet Cp—") 
sector and power law correlations in BTr/2, Xtt/2'^ 5) period-2 
VBS coexisting with one gapless mode in the singlet sector 
and power law correlations in 62fcFn '^2fcj?2' 6) phase with 
two gapless modes, one in the spin sector and one in the sin- 
glet sector. It is possible that some of the gapless phases will 
be further unstable to effects not considered here. 

The above essentially covers all natural possibilities of gap- 
ping out some or all of the low-energy modes of the generic 
SBM phase. Thus, as discussed at the end of Sec. IV A, we 
cannot pin (pa-+ because of the spin rotation invariance. The 
SU(2) spin invariance also imposes restrictions on the values 
of the variables that are pinned; these conditions are automati- 
cally satisfied in the above cases since our starting interactions 
are SU(2)-invariant. Furthermore, we cannot pin 9p^ because 
of the emergent conservation of J dx^Pp-- One exception is 
when the Fermi wavevectors take special commensurate val- 
ues; we discuss this next. 



D. Period-3 VBS state as a possible instability of the SBM in 
the commensurate case with kpi = n/3 

In the ring model Eq. (1), the DMRG observes translational 
symmetry breaking with period 3 in the intermediate parame- 
ter range flanked by the Spin Bose-Metal on both sides. Mo- 
tivated by this, we revisit the spinon-gauge theory in the spe- 
cial case with kpi ~ 7r/3 (then kp2 = — 57r/6). Compared 
to fermion interactions present for generic incommensurate 
Fermi wavevectors, we find one new allowed term 

Ve = U6(/^2T/fl2i/lla/i2T/L21./ma+H.C.) (60) 

= -Aue cos(V26'i^) sm{39p- - 9p+) . (61) 

The pinned 9p+ value is kept general at this stage. The scal- 
ing dimension is A[V6] = 1/2 + 9(7/4. Let us study what 
happens when g < 2/3 and Vg becomes relevant, so flows 
to large values. Then 6'p_ and 9icy are pinned while the conju- 
gate fields (/3p_ and (pi^ fluctuate wildly. There remains one 
gapless mode 6*20- that is still described by Eq. (20). 

We can use the bond energy and spin operators to character- 
ize the resulting state. First of all, B2kFi develops long-range 
order. Since 2kFi = 2tt/3, we thus have a Valence Bond 
Solid with period 3. Using Eqs. (A8) and (A14), the micro- 
scopic bond energy is 

6B{x) - cos(\/2e'i^) sin(2fcFix + kpi + 9p- + 9p+) .(62) 

In Eq. (61), we write 39 p^ - 9p+ = 3(6'p_ +9p+)- A9p+ and 
use the pinning condition on 46*^+, Eq. (17). There are two 
cases: 

a) uq cos{4:9 p-^) < 0: In this case, Vq is minimized by 
inequivalent pinning values V29ic ~ tt; 9p^ + 9p+ = 
7r/6,7r/6 + 27r/3,7r/6 + 47r/3. For 9p^ + 9p+ = 7r/6, the 
period-3 pattern of bonds is 

(56(a;) = {...,-l,i i ...}, (63) 

while the other two inequivalent pinning values give transla- 
tions of this pattern along the chain. A lower bond energy is 
interpreted as a stronger antiferromagnetic correlation on the 
bond. Then the above pattern corresponds to "dimerizing" ev- 
ery third bond as shown in Fig. 14a). 

b) U6cos(40p+) > 0: In this case, Vg is minimized by 
V29ia = 0; 9p^ + 9p+ = 7r/6,7r/6 + 27r/3,7r/6 + 47r/3. 
For 0p_ + 9p^ = 7r/6, the period-3 pattern of bonds is 

= -i ...}, (64) 

while the other two inequivalent pinning values give transla- 
tions of this along the chain. This pattern corresponds to every 
third bond being weaker as shown in Fig. 14b). 

Continuing with the characterization, we note that S2kFi 
and all operators at 7r/2 have exponentially decaying correla- 
tions. On the other hand, S2kF2 and B2kF2 have \/x power 
law correlations because of the remaining gapless 92c mode. 
Since 2fci?2 = 7r/3, we have period-6 spin correlations on the 
original ID chain. 
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a) 



kpi = 37r/8, we find a new quartic term 



b) 



FIG. 14: Valence Bond Solid states with period 3. Thick lines indi- 
cate stronger bonds. Remaining effective spin- 1/2 degrees of free- 
dom are also shown. Coexisting with the translational symmetry 
breaking, we have 1 /x power law spin correlations with the antifer- 
romagnetic (dynamic) pattern as shown. The two cases have slightly 
different microscopies, but are qualitatively similar on long length 
scales. 



The physical interpretation is simple. Consider first 
Fig. 14a) where every third bond is stronger A caricature 
of this state is that spins in the strong bonds form singlets 
and are effectively frozen out. The remaining "free" spins are 
separated by three lattice spacings and are weakly antiferro- 
magnetically coupled forming a new effective ID chain. Thus 
we naturally have Bethe-chain-like staggered spin and bond 
energy correlations in this subsystem, which coexist with the 
static period-3 VBS order in the whole system. The situa- 
tion in Fig. 14b) where every third bond is weaker is quali- 
tatively similar. Here we can associate an effective spin- 1/2 
with each three-site cluster formed by strong bonds. These 
effective spins are again separated by three lattice spacings 
and form a new weakly coupled Bethe chain. Note that while 
the theory analysis has the Fermi wavevectors tuned to the 
commensuration, the resulting state is a stable phase that can 
occupy a finite region in the parameter space, as found by the 
DMRG in the ring model (see Sec. V A). 

We can construct trial wavefunctions using spinons as fol- 
lows. In the mean field, we start with the band parameters ti 
and t2 such that kpi = tt/3 and then add period-3 modulation 
of the hoppings. The ±kpi Fermi points are connected by the 
modulation wavevector and are gapped out. The ±kF2 Fermi 
points remain gapless; just as in the Bethe chain case, the cor- 
responding Bosonized field theory provides an adequate de- 
scription of the long-wavelength physics, predicting 1/x de- 
cay of staggered spin and bond energy correlations. 

The above wavefunction construction and theoretical anal- 
ysis are implicitly in the regime where the residual spin cor- 
relations are antiferromagnetic. In a given physical system 
forming such a period-3 VBS, one can also imagine ferromag- 
netic residual interactions between the non-dimerized spins. 
Indeed, the DMRG finds some weak ferromagnetic tenden- 
cies in the ring model near the transition to this VBS state. 
This is not covered by our spin-singlet SBM theory but could 
possibly be covered starting with a partially polarized SBM 
state. 



E. Other possible commensurate points 

Alerted by the period-3 VBS case, we look for and find 
one additional commensurate case with an allowed new in- 
teraction that can destabilize the Spin Bose-Metal. When 



V4 = U4 .fl.^fl.^^^eal}fLla.fR2l3 + (R ^ L) + H.C. (65) 

~ -irii]V2] sin(26'p_ + 6'p+ - 9„+) sin((^p_ + Lp„^) 
-iViimi sin(26'p_ + dp+ + 6*^+) sin((y3p_ - ip^^) . 

(For the schematic writing here, we have ignored the commu- 
tations of the fields when separating the (p's and 9's.) The scal- 
ing dimension is A[V4] = 1/2 + g + l/(4g). This is smaller 
than 2 for g G (0.191, 1.309) and the interaction is relevant 
in this range. Since we have conjugate variables 6'p_ and (y9p_ 
both present in the above potential, we can not easily deter- 
mine the ultimate outcome of the runaway flow. It seems safe 
to assume that 6^+ and ipcr~ will be both pinned, which im- 
plies at least some period-2 translational symmetry breaking. 
One possibility, perhaps aided by the interactions Eq. (37), is 
that the ipp- is pinned; in this case, the situation is essentially 
the same as in Sec. IV B 2 and we get some period-4 struc- 
ture. Another possibility is that the 6'p_ is pinned; in this case 
{^ikpi) 7^ 0, and since Akpi = — 7r/2, we get period-4 bond 
pattern. 

To conclude, we note that the commensurate cases in this 
section and in Sec. IV D can be understood phenomenolog- 
ically by monitoring the wavevectors of the energy modes 
Bq. The dominant wavevectors are ±2fci?a, ±7r/2, and 
±4fci?2 = T4A:fi- When ^kp2 matches with 7r/2, we get 
the kpi = Stt/S commensuration of this section (here also 
2kpi matches with kp2 — kpi, while 2kp2 matches with 
3kpi + kp2). When Akp2 matches with 2kpi, we get the 
kpi = 7r/3 commensuration of the previous section. Tracking 
such singular wavevectors in the DMRG is then very helpful 
to alert us to possible commensuration instabilities, and both 
cases are realized in the ring model with additional antiferro- 
magnetic J3 = 0.5Ji discussed in Sec. VB, cf. Fig. 16. 



V. DMRG STUDY OF COMMENSURATE INSTABILITIES 
INSIDE THE SBM: VBS-3 AND CHIRALITY-4 

A. Valence Bond Solid with period 3 

As already mentioned in Sec. Ill, we find a range of pa- 
rameters where the SBM phase is unstable towards a Valence 
Bond Solid with a period of 3 lattice spacings (VBS-3). In the 
model with A'ling = Ji = 1, this occurs for 1.5 < J2 < 2.5, 
cf. Fig. 8. The characteristic correlations are shown in Fig. 15 
at a point J2 = 2. The dimer structure factor shows a Bragg 
peak at a wavevector 27r/3 corresponding to the period-3 VBS 
order. The spin structure factor has a singularity at a wavevec- 
tor 7r/3 corresponding to staggered correlations in the effec- 
tive spin- 1/2 chain formed by the non-dimerized spins, see 
Fig. 14. If we zoom in closer, the dimer structure factor also 
has a feature at 7r/3 that can be associated with this effective 
chain. 

To construct a trial VBS-3 wavefunction, we start with the 
spinon hopping problem that would produce kpi = 7r/3, so 
the first Fermi sea would be twice as large as the second. We 
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FIG. 15: (Color online) Spin and dimer structure factors at a rep- 
resentative point in the VBS-3 phase, J\ring = Ji = 1, J2 = 2, 
measured in the DMRG for system size L — 96 (we do not show 
the chirality as it is not informative.) The most notable features are 
the dimer Bragg peak at 2tt/3 corresponding to the static VBS order 
and also the spin singularity at tt/3 corresponding to the effective 
spin-1/2 chain formed by the non-dimerized spins, see Fig. 14. The 
trial VBS-3 wavefunction is constructed as described in the text. 



to 

of 

V 



o 

o jt/2 

CD 
> 



jr/3 



"ring 



: J.| = 1 , J3 = 0.5, varying J2; L=96 



m 
> 



a 

IE 
O 



Bethe, q^p 
SBM, q^igh 
SBM, qio„ 
VBS-3, q^F 
VBS-2 
„ Gutzw. 



♦ 



Belthe 




SBM 



^ VBS-3 



SBM VBS-2 



J2/J1 



FIG. 16: (Color online) Phase diagram of the Ji — J2 — A'ring 
model with additional antiferromagnetic third-neighbor coupling 
Ji — O.5J1 introduced to stabilize spin-singlet states. The study 
is along the same cut A'ring — Ji as in Fig. 8, and the overall fea- 
tures are similar, with the following differences. The ground state 
is singlet throughout eliminating the difficult "partial FM" region to 
the left of the VBS-3 (and the VBS-3 phase is somewhat wider). 
There is a sizable spin-gapped region between the Bethe-chain and 
SBM phases (see text for more details). A new spin-gapped phase 
with period-4 chirality order appears inside the SBM to the left of 
the VBS-3. 



in the dimer is stronger than the 2fci?i 2, which is in agreement 
with what we see in the neighboring SBM in Fig. 7. The re- 
emergence of the 4fci?2 and 2/sfi at the other end of the VBS- 
3 phase is obscured here by the weak ferromagnetic tendency 
(but is present in a model where this tendency is suppressed, 
see Fig. 16). 



Enhancement of the spin-singlet SBM by antiferromagnetic 
third-neighbor coupling J3 — 0.5 and a new phase with 
chirality order with period 4. 



then multiply every third first-neighbor hopping by 1 + (5 and 
Gutzwiller-project; for the point in Fig. 15 we find optimal 
S ~ 1. This gaps out the larger Fermi sea but leaves the 
smaller Fermi sea gapless. Our wavefunction is crude and 
shows a stronger dimer Bragg peak than the DMRG and some- 
what different spin correlations at short scales, but otherwise 
captures the qualitative features as can be seen in Fig. 15. 

The origin of the VBS-3 phase can be traced to the insta- 
bility of the SBM at special commensuration, cf. Sees. IV D- 
IV E. Indeed, in Fig. 7 we can follow the evolution of the 
singular wavevectors in the SBM phase between the VBS-2 
and VBS-3. As we decrease J2 moving towards the VBS- 
3, the 4fci?2 and 2kFi singular wavevectors in the bond en- 
ergy approach each other and coincide at 27r/3. When this 
happens, there is a new umklapp term that can destabilize the 
SBM and produce the VBS-3 state as analyzed in Sec. IV D. 
The instability requires g < 2/3 for the SBM Luttinger pa- 
rameter In this case according to Table I the 4fci?2 singularity 



As discussed in Sec. IIIC, in the original Ji — J2 — i^ring 
model, states in the SBM region near the left boundary of the 
VBS-3 tend to develop a small magnetic moment. We con- 
jecture that this occurs in the second spinon Fermi sea and 
suggest that an antiferromagnetic J3 will stabilize the SBM 
phase with spin-singlet ground state. One motivation comes 
from the picture of the neighboring VBS-3, where the non- 
dimerized spins are loosely associated with the second Fermi 
sea. These spins are three lattice spacings apart, so adding an- 
tiferromagnetic J3 should lead to stronger antiferromagnetic 
tendencies among them and also in the physics associated with 
the second Fermi sea. 

We have performed a detailed study adding a modest J3 = 
0.5 to the original model Eq. (1) along the same cut A'ling = 
Ji = 1. Our motivating expectations are indeed borne out. 
Figure 16 shows the phase diagram together with the evolution 
of the singular wavevectors as a function of J2/ Ji- While the 
overall features are similar to the phase diagram in the J3 = 
case. Fig. 8, a few points are worth mentioning. 
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First, the partial spin polarization is absent in the whole 
SBM phase between the Bethe-chain and VBS-3 phases. The 
DMRG converges confidently to spin-singlet ground state for 
L = 96. All properties are similar to those in Figs. 4 and 
6, providing further support for the singlet SBM phase in the 
original J3 = model. The VBS-3 phase and the SBM phase 
between the VBS-3 and VBS-2 are qualitatively very similar 
in the two cases J3 = and J3 = 0.5 and are not discussed 
further here. 

An interesting feature in the model with J3 = 0.5 is the 
presence of a sizable phase with spin gap intervening be- 
tween the Bethe-chain and SBM phases. Our best guess is 
that this phase has period-2 VBS order, although we do not 
see a clear signature in the dimer correlations. Our DMRG 
states in this region show very weak (if any) dimerization of 
the first-neighbor bonds, which may explain the lack of clear 
order in these dimer correlations. On the other hand, we see 
a sizable period-2 dimerization of the third-neighbor bonds 
{S{x) ■ S{x + 3)), but have not measured the corresponding 
bond-bond correlations to confirm long-range order We have 
not explored possible theoretical routes to understand the ori- 
gin of this phase yet and leave our discussion of this region as 
is. 

We now turn to one more new phase found in the J3 = 0.5 
model. In a narrow region inside the SBM phase not far 
from the left end of the VBS-3, we again find a spin-gapped 
phase. We identify this as having period-4 order in the chi- 
rality (Chirality-4 phase). Figure 17 presents a point J2 = 
1.5, J3 = 0.5. Looking at the singular wavevectors in Fig. 16, 
we see that this point is roughly where q\ow = 2k f2 passes 
IT /A. Analysis in Sec. IV E suggests an instability gapping out 
all modes and leading to some period-4 structure. Indeed, the 
DMRG spin and dimer structure factors show only some rem- 
nants of features near 2kFi, 2kp2, while the chirality shows 
a sharp peak at 7r/2. Looking at real space correlations, our 
L = 96 DMRG state breaks translational symmetry. The pat- 
tern of chirality correlations is consistent with the period-4 or- 
der shown in Fig. 13. The pattern of dimer correlations is also 
consistent with this picture and shows modulation with period 
2, which can be seen as a feature at tt in the dimer structure 
factor in Fig. 17. 

To summarize, with the help of modest J3 = 0.5 we have 
stabilized the spin-singlet SBM states between the Bethe- 
chain and VBS-3 phases. By suppressing potential weak fer- 
romagnetism, we have uncovered the Chirality-4 phase, which 
can be understood as arising from the instability of the SBM 
at the special commensuration discussed in Sec. IV E. 



VI. AN ATTEMPT TO BRING OUT PARTIALLY 
MAGNETIZED SPIN BOSE-METAL BY FERROMAGNETIC 
THIRD NEIGHBOR COUPLING J3 = 5 

In the original model Eq. (1), we do not have a clear un- 
derstanding of the "partial FM" states to the left of the VBS-3 
phase in Fig. 2. As discussed in Sec. IIIC, we suspect that 
there is a tendency to weak ferromagnetism in the second 
Fermi sea. However, the magnetizations that we measure are 
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FIG. 17: (Color online) Spin, dimer, and chirality structure factors 
at a point in the tentative Chirality-4 phase, A'ring = Ji = 1, J2 = 
1.5, J3 = 0.5, measured in the DMRG for system size L = 96. 
Note the absence of sharp features in the spin correlations, which 
suggests a spin gap, while the dimer correlations have only a feature 
at TV corresponding to period-2 modulation. On the other hand, the 
chirality structure factor shows a Bragg peak at 7r/2 that grows with 
increasing system size. The pattern of chirality correlations in real 
space is consistent with the order shown in Fig. 13. 
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not large: e.g., they are significantly smaller than if we were 
to fully polarize the second Fermi sea, and it is difficult to an- 
alyze such states. Here we seek better control over the spin 
polarization by adding modest ferromagnetic J3 = —0.5 in 
hopes of stabilizing states with a full spontaneous polariza- 
tion of the second Fermi sea, which is easier to analyze. 

We do not have as detailed phase diagram as for the J3 = 
and J3 = 0.5 cases. We expect it to look crudely similar to 
Figs. 8 and 16, with a narrower (if any) VBS-3 region, and 
with a wider partially polarized SBM region. We indeed find 
stronger ferromagnetic tendencies in the range < J2 < 2. 
However, we cannot claim achieving robust full polarization 
of the second Fermi sea and understanding all behavior The 
largest magnetization and properties closest to our expecta- 
tions are found in the middle region near J2 ~ 1. With this 
cautionary note warranting more work, we now present results 
for J2 ~ 1 to illustrate our thinking about such states. 

The DMRG study proceeds as follows. We start as before 
working in the S^^^^ = sector. We measure the spin structure 
factor and calculate the total spin ^tot from Eq. (32), which 
can give a first indication of a non-zero magnetization. How- 
ever, for the larger system sizes, the DMRG finds it difficult to 
converge to integer- valued S'tot due to a mixing of states with 
different total spins, and this leads to significant uncertainty. 
To check the value of the ground state spin, we run the DMRG 
in sectors with different 5^^^, expecting the ground state en- 
ergy to be the same for 5^^^ = 0, . . . , Stot and then to jump 
to a higher value for S^^^ > Stot- 

As an example, at a point J2 = 1 for system size L = 
48 we find that the DMRG energy is the same in the sectors 



C2 



0, . . . , 5 and then jumps, so the ground state spin is 



determined as St, 



5. The DMRG convergence is good 



and the SU(2)-invariant structure factor {Sq ■ S-q) is the same 
measured in the different sectors S^^^ < 5, indicating that 
these states indeed belong to the same multiplet. 

The situation is less clear for i = 96 because of reduced 
convergence. At the point J2 = 1 the extensive energies 
obtained by the DMRG in the sectors Stot = 0, . . . , 10 are 
non-systematic and are within O.lJi of each other, with the 
lowest energy found in the S^^t — 10 sector. Also, the SU(2)- 
invariant structure factors differ slightly and the estimates of 
5tot vary around S'tot ~ 8 — 10. The convergence is best in 
the highest S^^^ = 10 sector, where the total spin is found to 
be accurately S'tot = 10; the improved convergence is indeed 
expected since there are fewer available low energy excited 
states to mix with (e.g., spin-wave excitations of the ferro- 
magnet are not present in the highest S^ sector). Interestingly, 
we find a fully converged state in the S'tot = 11 sector with 
Stot = 11 whose energy is only slightly higher, which prob- 
ably adds to the above convergence difficulties. More impor- 
tantly, the energy jumps to a significantly higher value in the 
sector Sfot = 12- Our best conclusion is that the total spin of 
the ground state is Stot = 10. 

Turning to the VMC study, we consider a family of varia- 
tional Gutzwiller states where we allow different spin up and 
spin down populations of the two Fermi seas centered around 
k = and tt (we do not attempt any further improvements on 
top of such bare wavef unctions). In the model parameter re- 
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FIG. 18: (Color online) Spin and dimer structure factors at a tenta- 
tive point with partial ferromagnetism, A'ring — Ji = 1, J2 = 1, 
J3 = —0.5, measured in the DMRG for system size L = 96. The 
calculations are done in the sector Stot = 10 where the DMRG is 
well-converged and gives S'tot ~ 10, which we think is the true 
ground state spin. The VMC state has the second Fermi sea fully 
polarized with N2-\ = 20, N2i = 0, while the first Fermi sea is 
unpolarized with Ni^- — Nii — 38. Vertical lines label important 
wavevectors 2kFi, 2kF2, and —kF2 — kpi- 



gion discussed here, we find that the optimal such states have 
a fully polarized second Fermi sea and an unpolarized first 
Fermi sea. For the L ~ 48 example quoted above, the optimal 
VMC polarization indeed matches the DMRG Stot = 5, while 
for the L = 96 case the optimal VMC state has Stot = 11 and 
a state with Stot = 10 is very close in energy. Appendix B 3 
provides more details on the properties of such Gutzwiller 
states, while here we simply compare the VMC and DMRG 
measurements. 

The DMRG structure factors for the J2 = 1, J3 = -0.5, 
L = 96 system are shown in Fig. 18, together with the VMC 
results for the Gutzwiller wavefunction with 5tot = 10. A no- 
table difference from the singlet SBM states of Sec. Ill is that 
the characteristic peaks are no longer located symmetrically 
about tt/2. For example, in the trial state we have prominent 
wavevectors 2kp2 and 2kFi that satisfy 2kp2 + 2 x 2kFi = 
2tt; also, we have a wavevector —kp2 — kpi, which is now 
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FIG. 19: (Color online) Separate (SgS'^g) and (SgSlg) structure 
factors for the same system as in Fig. 18. 



different from tt/2. It is not easy to discern all wavevectors 
in Fig. 18 because the 2fci?2 happens to be close with the 
—kp2 — kpi. Nevertheless, the overall match between the 
DMRG and VMC suggests that the trial wavefunction cap- 
tures reasonably the nature of the ground state. 

Working in the sector S'^^j ~ Stot also allows more de- 
tailed comparison between the DMRG and VMC. In this case, 
there is a sharp distinction between the (S^Stq) and (S'^S'f ^) 
structure factors. The former has singular wavevectors 2kpi, 
2kF2, and —kp2—kpi, while the latter is lacking the wavevec- 
tor 2kp2 since there is no spin-flip process across the second 
Fermi sea. Our measurements are shown in Fig. 19. The 
wavevectors 2k p2 and — fcF2 — kpi are too close to make a 
more clear-cut distinction; nevertheless, the VMC reproduces 
all details quite well. 

In analogy with the SBM theory, we expect the 2k pi and 
2k p2 singularities to become stronger compared with the bare 
Gutzwiller and the —kp2 — kpi to become weaker; this is 
roughly consistent with what we see in the DMRG structure 
factors. As discussed in Appendix B 3, however, we do not 
have a complete description of such a partially polarized SBM 
phase that must incorporate ferromagnetic spin waves as well 
as the low energy SBM modes. This is left for future work. 



We also mention that the closeness of the 2kp2 and ~kp2 — 
kpi warns us that the system is near a commensuration point 
with kpi = 2tt/5 where it can be further unstable, which 
requires more study. 

To summarize, by adding modest ferromagnetic J3 = —0.5 
we have reaUzed the SBM state with fully polarized sec- 
ond Fermi sea, confirming our intuition about the origin of 
the weak ferromagnetic tendencies in the original model dis- 
cussed in Sec. Ill C. A more thorough exploration of the phase 
diagram in the model with J3 = —0.5 as well as in the original 
model in the partial FM region is clearly warranted to develop 
better understanding of such partially ferromagnetic states. 



VII. CONCLUSIONS AND FUTURE DIRECTIONS 

We have summarized the main results and presented much 
discussion particularly in Sec. I. Perhaps one point we would 
like to reiterate is the remarkable coincidence between the 
sign structure present in the DMRG wavefunctions for the 
spin model SBM phase and the sign structure in the spin sector 
of free fermions on the ladder (e.g., metallic electrons). This 
sign structure is encoded in the singular wavevectors ("Bose 
surfaces"), and indeed the Gutzwiller-projected wavefunction 
with just one variational parameter is sufficient to reproduce 
the locations of all the singularities throughout the observed 
SBM phase. 

We conclude by mentioning some standing questions and 
future directions. First, in the ring model, we have focused 
on the Spin Bose-Metal and dealt with other phases only as 
needed to sketch the rich phase diagram. For example, we 
have not studied carefully the VBS-2 region, which might 
harbor additional phases. We have not studied adequately 
(numerically or analytically) the spin-gap region between the 
Bethe-chain and SBM phases in the J3 ~ 0.5 model. Fig. 16. 
One question to ask here is whether there is a generic instabil- 
ity when we start populating the second Fermi sea, or whether 
we can go directly from the Bethe-chain phase to the SBM. 
More generally, we have not studied various phase transitions 
in the system. 

Next, while we understand the long-wavelength SBM the- 
ory with its single Luttinger parameter g, we have found that 
the Gutzwiller wavefunctions represent only the special case 
g = 1 and cannot capture the general situation g < 1. It 
would clearly be desirable to construct spin-singlet wavefunc- 
tions appropriate for the general case. Even thinking about 
the Gutzwiller wavefunctions, it could be interesting to un- 
derstand the observed g — I analytically and ask if they may 
be exact ground states of some Hamiltonians (in the spirit of 
the Haldane-Shastry model^^'^^). 

On a separate front, we have encountered an interesting 
possibility of the Spin Bose-Metal with partial ferromag- 
netism occurring in one of the subbands, but more work is 
needed to fully understand the numerical observations and de- 
velop analytical theory. Even without a spontaneous moment 
in the ground state, our observations suggest that in the regime 
between the Bethe-chain and VBS-3 phases the second band 
is narrow in energy. Some ferromagnetic instability or possi- 
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bility of spin-incoherent regime in one of the bands can lead 
to anomalous transport properties in such a quantum wire, a 
topic of much current interest. 

Looking into future, it will be interesting to consider elec- 
tron Hubbard-like models on the two-leg triangular strip and 
look for possible SBM phase. The Hubbard model has been 
studied in a number of works,'*^"^'^ but the focus has been 
mainly on the conventional insulating phases such as the 
Bethe-chain and VBS-2 states. This is appropriate in the 
strong Mott insulator limit, ii,t2 <C U, where the effec- 
tive spin model is the Ji — J2 model with ,/i = 2t\/U , 
J2 = 2t^/U . However, at intermediate coupling just on 
the insulator side, one needs to include multiple-spin ex- 
changes, and the leading new term is the ring exchange with 
i^ring = 2{)t\t\/U^ . As we havc learned, this ring term sta- 
bilizes the SBM phase, so revisiting the Hubbard model with 
the insights gained here is promising. In Sec. II B, we ap- 
proached the SBM by starting with a metallic two-band elec- 
tron system ("C2S2") and gapping out only the overall charge 
mode Then it is natural to look for the SBM near an 

extended such C2S2 metallic phase, and we may need to con- 
sider electron models with further-neighbor repulsion to open 
wider windows of such phases. 

Last but not least, we would like to advance the program 
of ladder studies closer to 2D. It is prudent to focus on the 
spin model with ring exchanges. On the exact numerics front, 
4 to 6 legs is probably at the limit of the DMRG capabili- 
ties. The VMC approach should still be able to capture the 
critical surfaces if they are present, since they are dictated by 
short-distance physics; on the other hand, the bare Gutzwiller 
will likely fail even more in reproducing correct long-distance 
behavior. We do not know how far the present Bosonization 
approach can sensibly hold going to more legs. These are 
challenging but worthwhile endeavors given the experimental 
importance of understanding weak Mott insulators. 
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APPENDIX A: OBSERVABLES IN THE SBM PHASE 

We have defined spin S(x), bond energy B{x), and chiral- 
ity x{^) observables in Sec. II C [cf. Eqs. (22-23)]. Here we 
find detailed Bosonized forms by a systematic construction of 
observables in the SBM (and will find more observables on 
the way). 

In the gauge theory treatment Sec. II A, we consider gauge- 
invariant objects constructed from the spinon fields; in the in- 
teracting electron picture of Sec. II B, these are operators that 



do not change the total charge. 

We begin with fermion bilinears and first consider the ones 
composed of a particle and a hole moving in opposite di- 
rections. Such bilinears are expected to be enhanced by 
gauge fluctuations since parallel gauge currents experience 
Amperean attraction. We organize these bilinears as follows: 

S2kF^ = -^flaa^apfRalS , (Al) 

^2kFa = 2-^LaafRaa , (A2) 

Sn/2 = ■^fRia^afifL2P + ^fl2a^aPfLll3 , (A3) 

<^7t/2 = ■^fRiafL2a + -j^fmjLla , (A4) 

4/2 = ■^fRia^al3fL2P - -^fR2a^apfLip , (A5) 

X7T/2 = -^fRiafL2a - -^fR2cjLla , (A6) 

with 5_Q = Sq, etc. The microscopic spin operator ex- 
panded in terms of the continuum fermion fields readily gives 
the listed Sq. 

The bond energy can be approximated as the spinon hop- 
ping energy, 

B{x)^-t[fi{x)Ux + l) + R.c.] (A7) 

(recall that we work in the gauge with zero spatial vector po- 
tential). Expansion in terms of the continuum fields gives, up 
to a real factor, 

Bq - e'^^/^eg . (A8) 

Such connection between Bq and eq is understood below 
for all Q ^ TT. The objects eg are convenient because of 
their simpler transformation properties under lattice inversion. 

The physical meaning of the operators 5^/2 ™d Xtt/2 can be 
established on symmetry grounds. Thus, x is the spin chirality 
defined in Eq. (23). The expression in terms of bilinears can 
be also found directly by considering the circulation of the 
gauge charge current around the [x — l, x, x + 1] loop, 

Xix) -^itrAflir)Ur') - H.c] . (A9) 

This is familiar in slave particle treatments: the circulation 
produces internal gauge flux whose physical meaning is the 
spin chirality. On the other hand, S is related to the following 
microscopic operator; 

D{x) = S{x) X Six + 1) , D^/2 = e"/''4/2 ■ (AlO) 

At Q — tt/2, this enters on par with S, B, and x- 
The bosonized expressions at the 2kFa are: 

S2k,^ - -iVa^Vaie^'^+e^^'^- sm{V2cpaa) , (All) 

S^kF. ^ -^Va^Vaie'^'+e^'''^- cos{V2^aa) , (A12) 

^2fe,„ = -e*«''+e±'^-sin(V20„,), (A13) 

e2fc^„ = ie^'^+e^'O^-cosiv^eaa) , (A14) 
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where the upper or lower sign in the exponent corresponds to 
a = 1 or 2. The pinned value 9p^, which is determined by 
minimizing Eq. (17), is left general at this stage. It is not so 
important for the qualitative behavior at the 2kpa and tt/2, 
but is crucial at a wavevector tt later. 

For each a, the {S, B)2kFa structure is similar to that in a 
single Bethe chain except for the 9p exponentials. As in the 
Bethe chain, we expect the spin and VBS correlations to be 
closely related - in particular, they decay with the same power 
law. The corresponding scaling dimension in the fixed-point 
theory Eq. (18) is 



A[S2fc.J = i + f 



The bosonized expressions at the tt/2 are: 



(A15) 



X7T/2 



i77i|?72xe 



sin(v3p_ + ip„+) (A16) 



iViim^e''^''' sin((^p_ - ip^+) 



(A17) 



, (A19) 
VJa_)(A20) 



(A21) 



sin((y9p_ + ipa-) (All) 



(A23) 



Vi]V2te 



C0s(v3p_ + (y3^_) (A24) 



(A25) 



Expressions for 8.,^ 12 can be obtained from those for es- 
sentially by interchanging sines and cosines. As before, Bt^i2 
is given by Eq. (AS). The above details are needed partic- 
ularly when we discuss phases arising as instabilities of the 
SBM, Sees. IV B-IV E, while in the SBM we immediately see 
that all scaling dimensions are equal: 

(A26) 

This completes the "enhanced" bilinears. We also mention, 
without giving detailed expressions, "non-enhanced" bilinears 
at wavevectors Q = ±(fcF2 — kpi)- Their scaling dimension 
is 

A[Sq] - A[Bq] = A[Dq] = A[xq] = ^ + ^ + f , (A27) 

which is always larger than the spinon mean field value of 1. 

Finally, we have bilinears carrying zero momentum - es- 
sentially JpaaiJpaa from Eq. (34). These give conserved 
densities and currents and have scaling dimension 1. We 



specifically mention examples leading to Eqs. (25-27): 

Sq=0 ^ Jrii + Jill + '^fl22 + "^£22 — —dx&a+ , (A28) 

TT 

2 

£0=0 ^ Jrii + Jlii — Jb.22 — Jl22 = —dxOp- , (A29) 

TT 

2 

Xq=o ~ Jbii - Jlii - Jr22 + Jl22 = -d^^Pp^ .(A30) 

TT 

(One way we can make the identifications in the last two lines 
is by using physical symmetry arguments.) 

So far, we have only considered fermion bilinears. Since 
the theory is strongly coupled, we should also study contribu- 
tions with more fermion fields. We now include four-fermion 
terms focusing on the spin, bond energy, and chirality oper- 
ators that are measured in the DMRG. First, there appears a 
new wavevector Akpi = —AkF2 in the bond energy, via, 

e4fc,, : flnflnfRnfmi - e'^'^+e^''^- , (A31) 

/K2T/«2i/i2T/L2i ^ e-'2e,+ g.2e,_ (^32) 

The two contributions come with independent numerical fac- 
tors and can be also generated as {e2kFi)'^ ^"d (e_2fcp2)^- 
Once the 9p+ is pinned, there is only one qualitatively distinct 
contribution and the scaling dimension is 



A[BikFr]=9 



(A33) 



Note that for sufficiently small g < 2 /3, the power law decay 
is slower than that of the bilinears B2kFa ■ There is no compa- 
rable Akpi contribution to the spin operator. 

Four-fermion terms bring out another important wavevec- 
tor, Q = TT. We list independent dominant such contributions 
to S'tt, Btt, and Xn- 

5^: sin(20,+ )sin(20p+) , sin(20<,„) sin(20p+) ;(A34) 
B^ : [008(26*^+) + cos(26l^_)] sin(26lp+) , (A35) 
[cos(26'<,+ ) + f cos(2(/j^_)] sin(26lp+) , (A36) 
fcos(2(pp_)sin(26'p+) ; (A37) 
X^: fsin(2(/jp_)sin(20p+) . (A38) 

These can be generated by combining the previously exhibited 
bilinears as follows. 5*^: S2i^^_^e2kF2^^2kF2^'^'^f"^ 

TT- 

ie2kFi^2kF2 + H.C., i(e2^2 - xl/2) + H.c, ^(e^^a + xl/2) + 
H.c; Xtt: X-n/2^-n/2 + H.c. The scaling dimensions are 



A[S^] = A[B^] = 1 
A[x.] = 1/5. 



(A39) 
(A40) 



The above observables are present if sin(26'p_|-) 7^ 0, e.g., if 
the Op+ is pinned as in Eq. (14), which we argued is natu- 
ral when the spin model is describing a Mott insulator phase 
of a repulsive electron model. On the other hand, the above 
contributions would vanish if the pinning potential Eq. (17) 
had i'8 < 0. Some other physical observables containing 
cos(26'p+) and having different symmetry properties would 
be present instead. We do not write these out since both the 
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DMRG and the trial wavefunctions have signatures in the spin, 
VBS, and chirality at the wavevector tt, suggesting that the 
pinning Eq. (14) is realized. We have identified several more 
new observables at tt containing sm(29p+); we do not spell 
these out here since our primary focus is to understand fea- 
tures in the numerics measuring the familiar S, B, and x- 

We finally mention that four-fermion terms produce still 
more wavevectors, ±(3/o>i + kp2) = T(3fcF2 + ^Fi); for ex- 
ample, ikpi + kF2 can be obtained by combining 2kFi and 
— 7r/2. The scaling dimensions are the same as at ±(A;f2 — 
A:Fi),Eq. (A27). 

For completeness, we have also checked six-fermion and 
eight-fermion terms. The only new wavevectors where the 
scaling dimension can be smaller than 2 are QUpa = '^kpa + 
^kpa (scaling dimension 1/2 + 9.g/4), and Skpa = ^kpa + 
Akpa (scaling dimension Ag). However, one needs small g for 
these to become visible and in any case they always have faster 
power law decay than at 2fcFa and Akpa- Finally, entries listed 
as "subd." in Table I can be constructed, e.g., as XQ ~ XocQj 
which has scaling dimensions 1 + A [eg]. 

Table I summarizes our results for the correlations in the 
Spin Bose-Metal phase. In words, we expect dominant spin 
and VBS correlations at the wavevectors ±2fcFi, ±2fcF2 de- 
caying as l/x^+^Z^ and at the wavevectors ±7r/2 decaying as 
l/x^+^Z*^^^'. The former decay is more slow since stability 
of the phase requires g < \. Note that the wavevectors 2fcFi 
and 2kF2 are located symmetrically around tt/2. We also ex- 
pect a power law at the wavevectors and tt. Next, at 
the wavevectors ±(A:f2 — kpi) and ±(3fcFi + kpT), which 
are also located symmetrically around ±7r/2, we expect a still 
faster power law l/x^+^/'^^^'+^Z^. Furthermore, the bond en- 
ergy shows a power law X/x^^ at -izAkpi. The spin chirality 
has similar signatures to the above at ±7r/2, 0, ±(fcF2 ~ ^Fi), 
and ±(3/cfi + kF2), but decays as Xjx^l^ at tt. These are 
the simplest observables that can be used to identify the SBM 
phase in a given system. Figure 21 shows measurements in the 
Gutzwiller wavefunction projecting two Fermi seas and nicely 
illustrates all singular wavevectors, while it appears that such 
wavefunctions reahze a special case with 5 = 1. We also re- 
mark that in the general SBM the presence of the marginally 
irrelevant interactions Eq. (43) will lead to logarithmic correc- 
tions in correlations.^''^- 

We conclude by describing our treatment of the Klein fac- 
tors (see, e.g., Ref. 32 for more details). We need this when 
determining "order parameters" of various phases obtained as 
instabihties of the SBM, Sees. IVB-IVE. The operator f 
from Eq. (40) has eigenvalues ±1. For concreteness, we work 
with the eigenstate corresponding to +1: r|+) = |+). We 
then find the following relation 



|77l-|-772t| 



l'7u'72il 



pure imaginary . (A41) 



wavevectors; for example, 



£^2 -e ''+(+I'UT'/2tI 



cos((^p_) sin(6'CT+) sin((/3CT_) 



+ i sin(v3p_ ) cos(6'ct+) cos((/3cr-) 



, (A42) 



X7r/2 = -e '^'■+(+|?7it?72tI+) cos((y9p_)cos(6l^+)cos(v?^_) 



+ I sin(v3p_ ) sin(6'a+ )sm{^Pa~) 



(A43) 



APPENDIX B: DETAILS OF THE WAVEFUNCTIONS 
1. Gutzwiller projection of two Fermi seas 

It is convenient to view the spin wavefunction as that 
of hard-core bosons, where up/down spin corresponds to 
present/absent boson. In the general spinon construction, 
we occupy = 1, • • • , ^ti orbitals with spin up and 

{/cj, j = 1, . . . , TVj^} orbitals with spin down; iV| + -V^ = L 
is the size of the system. After the Gutzwiller projection, the 
boson wavefunction is 



z = 1, . . . , 7V^}) = detie'^^l^'] det[e^^'^«-] ,(B1) 



where the set {pj , j = 1 , . . . , A^i } is a complement to 



{_fcj, j 1, . . . , iVx} in the Brilloin Zone (BZ). The mo- 
mentum carried by this wavefunction is 'Y^=\^j + Pi) = 

TV T N I 

Ej=i + Ejii + EgGBz 9- particular, we see that 
the wavefunction remains unchanged if we shift all occupied 
spinon momenta by the same integer multiple of 2tt/L. 



We now consider the spin-singlet case when 



N 



This is useful when discussing observables at the ±7r/2 



L/2, {k^} = {fcJ-}. Here L is even and all momenta are in- 
teger multiples of 27r/L, so Egesz 9 = For convenience, 
we assume that L is a multiple of 4. Fig. 20 illustrates filled 
fc-points for the band in Fig. 3. We have two Fermi seas of vol- 
ume A^i and N2 in the symmetric configuration, i.e., separated 
by L/4 unoccupied orbitals on each side around the Brilloin 
Zone. Relating to the band Fig. 3, the larger iVi corresponds 
to occupied /c -points centered around 0, while N2 corresponds 
to points around tt. As already noted, a solid shift of the occu- 
pied states leaves the wavefunction unchanged. We can then 
specify such symmetric state as (iVi, A^2), which requires only 
one parameter since Ni + N2 = = N ^ = L/2. We can 
readily verify that such {even, even) states carry momentum 
and are even under site inversion operation while {odd, odd) 
states carry momentum tt and are odd under inversion. 

The relative wavevectors connecting the Fermi points are 
gauge-independent and are observed in various quantities, see 
Fig. 21. Specifically, using variational Monte Carlo, ^'^•"'" we 
measure the spin structure factor and see dominant singulari- 
ties at the wavevectors ±2fcFi, ±2kp2, and ±{kpi + kp2) = 
=F7r/2, which connect Fermi points with opposite group ve- 
locities. By studying sizes up to L = 512 and performing 
scaling analysis at these wavevectors, the singularities appear 
to have the same power law. This is consistent only with the 
special case g — 1 in the SBM theory [cf. Eqs. (A15, A26) 
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Bare Gutzwiller: L=256, (^^,^2) = (104,24) 



FIG. 20: View of the wavefunction constructed by filling fc-states of 
spinons. Here momenta k = 2nn/L, n = 0, 1, . . . , L — 1, form a 
closed circle. Each filled dot is occupied by both spin up and spin 
down, producing spin singlet. The projected wavefunction remains 
unchanged if all momenta are shifted by the same amount. Only the 
relative configuration matters, and here we show symmetric config- 
uration of the two Fermi seas separated by L/4 unoccupied fc-states 
on either side. This is our "bare Gutzwiller" wavefunction for the 
Spin Bose-Metal. 



and Table I]. Our direct estimates of the scaling dimensions 
are also consistent with the value A = 3/4 expected in this 
case. Turning to other less singular wavevectors, we clearly 
see V-shaped (~ \5q\) features at and tt corresponding to 
scaling dimension 1, which is expected generally. We also 
see ±(fci?2 — kpi) with scaling dimension 1, which requires 

We next consider VBS correlations and see all of the above 
wavevectors, but we can not quantify the singularities as ac- 
curately. The VBS correlations also show singularities at 
±(3fci?i + fci?2) and ±4fci?2 = T4fcFi (the former is also 
expected in the spin structure factor but is not visible there, 
probably due to amplitude effect, while the ±4fci?2 is expected 
only in the bond energy). 

Finally, we measure spin chirality correlations and see 
dominant singularity at ±7r/2 consistent with A — 3/4. We 
also see singularities at wavevectors 0, tt, ±{kF2 — kpi), and 
±(3fcFi + kF2) consistent with A = 1, again as expected in 
the special case g = I- 

The above appears to hold for a range of relative popula- 
tions of the Fermi seas (away from the limiting situations of 
a single or two equal Fermi seas). We are then led to conjec- 
ture that such spin-singlet wavefunctions with two Fermi seas 
have correlations given by the SBM theory of Sec. II and Ap- 
pendix A with .9 = 1. This conjecture is natural since in the 
theory the parameter g depends on the ratio of the two Fermi 
velocities, cf. Eq. (21), while the wavefunction knows only 
about the occupied/unoccupied states and does not contain the 
band energy parameters. We leave proving this conjecture an- 
alytically as an open problem. 

Given the preceding discussion, it appears that such bare 
Gutzwiller wavefunctions cannot capture fully the properties 
of the generic Spin Bose-Metal as described by the theory of 
Sec. II with general g < 1. It is possible that they are appro- 
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FIG. 21 : Spin, bond energy, and chirality structure factors in the bare 
Gutzwiller wavefunction with two Fermi seas (A'^i , N2) — (104, 24) 
on the ID chain of length L = 256. The expected singular wavevec- 
tors are marked by vertical lines (here kpi and fcF2 are defined as in 
Fig. 3 and all indicated wavevectors are modulo 27r). The structure 
factors are symmetric with respect to g —q, and we only show 
< 5 < TT. The character of the singularities is consistent with the 
special case <; = 1 in the SBM theory of Sec. II. 
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priate wavefunctions for some critical end-points of the SBM 
phase where the parameter g = 1, e.g., for the transitions out 
of the SBM discussed in Sec. IV B. While the energetics study 
with the bare Gutzwiller wavefunctions gives us first indica- 
tions for the SBM phase, it is desirable to have more accurate 
trial states. This is what we turn to next, although only with 
limited success. 



2. SU(2)-invariant improvement of the Gutzwiller 
wavefunctions 

To allow more variational freedom, we consider mean field 
with both spinon hopping ^(fc), Eq. (3), and spinon pairing in 
the singlet channel with real gap function A{k) (this way, the 
wavefunction remains spin rotation and time reversal invari- 
ant). Generic A(fc) would open up gaps, while we want the 
wavefunction to be critical. One way to maintain gaplessness 
is to require A(fc) to vanish at the Fermi points. This can be 
achieved, for example, by taking 

A(fc) = fimk) (B2) 

with a smooth /(fc). We have tried several simple functions 
f{k), e.g., expanding in harmonics 

= I]/«cos(nfc) , (B3) 

n 

with few /„ treated as variational parameters. Upon writing 
out the corresponding Gutzwiller wavefunction, one can see 
that the dispersion ^(fc) enters only through its sign ^ < or 
^ > 0, so in the case with two Fermi seas like in Fig. 20 we 
can use the same label {Ni, N2) and expect similar singular 
wavevectors encoded in the relative positions of the "Fermi 
points." We refer to such a state as "improved Gutzwiller" 
and can view it as a "gapless superconductor," although with 
caution because of the non-intuitive effects of the projection. 
For example, the SU(2) gauge structure of the projective con- 
struction implies that /„ — A5n,o gives the same state as the 
bare Gutzwiller independent of A. In practice, we often fix /q 
and vary /i,/2. 

For the zigzag ring model in the Spin Bose-Metal regime, 
such approach improves the trial energy by about 50-60% 
compared to the difference between the bare Gutzwiller en- 
ergy and the exact DMRG ground state energy. The exponents 
of the power law correlations in the improved wavefunctions 
appear to remain unchanged from the bare case, although the 
numerical amplitudes are redistributed to resemble the DMRG 
correlations better, as can be seen in the examples in Sec. Ill B, 
Figs. 4 and 5. Thus, this approach is only partially success- 
ful since we can not produce the long-distance behavior ex- 
pected in the generic SBM and tentatively seen in the DMRG. 
Still, the fact that we can significantly improve the trial energy 
while retaining the underlying gapless character gives us more 
confidence in the variational identification of the SBM phase. 

We also mention that in the Bethe-chain regime where the 
bare Gutzwiller projects one Fermi sea, the "gapless super- 
conductor" improvement with parameters /g , /2 works even 



better, bringing the trial energy much closer to the exact 
DMRG value and better reproducing short-scale features in 
the spin correlations, see Fig. 22. Such good trial states for 
the competing Bethe-chain phase give our VMC more accu- 
racy in determining where the SBM phase wins energetically 
and more confidence interpreting the DMRG results. 



3. States with fully polarized second Fermi sea 

Motivated by the possibility of partial ferromagnetism in 
some regimes discovered in the DMRG study of the ring 
model (see in particular Sec. VI), we have also considered 
Gutzwiller projection of states with unpolarized large Fermi 
sea and fully polarized small Fermi sea. The spin correla- 
tions here can be understood using a naive Bosonization treat- 
ment starting with such spinon mean field state and following 
the same procedure as for the unpolarized Spin Bose-Metal 
in Sec. II. The naive long-wavelength theory now has two 
free Boson modes. The dominant correlations are expected 
to be at wavevectors that connect Fermi points with oppo- 
site group velocities. Taking the polarization axis to be z, 
the spin structure factor {SgStq) has dominant singularities 
at 2fcFi, 2fci?2, and —kF2 — kpi, while the {S^S^^) is miss- 
ing the 2kp2 since there is no spin-flip process across the sec- 
ond (polarized) Fermi sea. We indeed observe such correla- 
tions in the wavefunctions, and the dominant power law enve- 
lope is consistent with x^^^"^, which is what such naive theory 
would give if we assume equal velocities near all Fermi points 
and ignore all interactions other than gapping out the overall 
"charge" mode, Eq. (13). We note, however, that to prop- 
erly describe such a partially polarized phase in the system 
with short-range interactions, we would need to also account 
for the ferromagnetic spin wave, which is not present in our 
wavefunctions"*- and not treated in the more general (but still 
naive) Bosonization theory outlined above. We do not pursue 
this further here. 



APPENDIX C: DMRG RESULTS IN CONVENTIONAL 
PHASES ON THE ZIGZAG CHAIN 

For ease of comparisons, here we show our DMRG mea- 
surements in the conventional Bethe-chain and VBS-2 phases 
identified on the zigzag chain in earlier works. We take 
representative points from the same cut ifiing/'/i = 1 studied 
in detail in Sec. Ill, since this allows us to better relate to the 
SBM phase at such significant A'ling values. 



1. Bethe-chain phase 

Figure 22 shows spin, dimer, and chirality structure fac- 
tors at J2 = —1, measured in the DMRG for system size 
L = 192. The DMRG can still obtain reliable results with 
m = 3200 states kept in each block, and this is related to 
the smaller central charge than in the SBM phase (as dis- 
cussed in Sec. HID). Figure 22 also shows the structure fac- 
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FIG. 22: (Color online) Spin, dimer, and chirality structure factors 
at a representative point in the Bethe-chain phase, A'ring = Ji — 1, 
J2 — —1, measured in the DMRG for system size L = 192. We 
also show structure factors in the bare Gutzwiller-projected single 
Fermi sea state, {Ni,N2) = (96, 0), and in the improved Gutzwiller 
wavefunction with parameters /o = l,/i = 0, /2 = —1.4 (see 
Appendix B for wavefunction details). 



tors in the bare Gutzwiller-projected single Fermi sea state, 
{Ni,N2) = (96,0), and in the improved Gutzwiller wave- 
function (see Appendix B); the latter achieves significantly 
better trial energy and overall match with the DMRG results. 

At long distances, we expect both the spin and dimer cor- 
relations to decay with the same power law: {S{x) ■ S{Q)) ^ 
{—ly^/x, {B{x)B{Q)) ^ (—l)^7a;, up to logarithmic correc- 
tions. We indeed see roughly such power law in the real space 
correlations. Some quantitative aspects are different from the 
pure Heisenberg chain, for which the bare Gutzwiller state 
is a good approximation. Thus, the spin structure factor in 
Fig. 22 has a larger amplitude of the g = tt singularity and 
also develops a hump at wavevectors below 7r/2. Both these 
features are captured by the improved Gutzwiller wavefunc- 
tion. On the other hand, the dimer structure factor has a sig- 
nificantly smaller amplitude of the q = tt singularity than 
the pure Heisenberg chain and the bare Gutzwiller; the im- 
proved Gutzwiller wavefunction moves in the right direction 
compared to the bare one but still does not capture well the 
amplitude at tt. 

For the chirality correlations, we expect {x{x)x{0)) ^ 
{—l)^/x^ + l/x'*, and we indeed see some fast decay in the 
real space data comparable with the power law behavior. The 
corresponding momentum space singularities at g = tt and 
g = are very weak. In agreement with this, we do not see 
any features in the chirality structure factor in Fig. 22. 

This Bethe-chain phase example allows to contrast with 
the SBM phase in Sec. IIIB, where we see different singu- 
lar wavevectors and prominent features in all these observ- 
ables including the chirality. The experience of being able to 
improve significantly the short-scale features in the trial wave- 
functions carries over to the SBM, although in the Bethe-chain 
phase we have an advantage that our wavefunctions also cap- 
ture the long-distance power laws correctly. 



2. Valence Bond Solid with period 2 

Consider now the large J2 case. In the J2 ^ 00 limit, we 
have decoupled legs, and each behaves as a Heisenberg spin 
chain. Finite J1/J2 and K^i^g/ J2 will couple the two legs and 
will likely open a spin gap'^ '^ producing a VBS state with pe- 
riod 2 (Fig. 11). Figure 23 shows our measurements at a repre- 
sentative point J2 = 4 from the Kiing = Ji = 1 cut. The spin 
correlations show a dominant peak at a wavevector q ^ tt/2 
and bond correlations have a peak at q = it. We compare with 
the Gutzwiller projection of two equal Fermi seas in the ID 
zigzag chain language, or, equivalently, decoupled legs in the 
two-leg ladder picture. This wavefunction is thus strictly ap- 
propriate only in the J2 — > 00 limit, but it clearly reproduces 
the DMRG data quite well. 

Looking at Fig. 23, there is not much direct evidence for 
the VBS-2 order in the DMRG data. It is safe to say only that 
upon exiting the SBM phase along this cut, we are close to 
the fixed point of decoupled legs. One argument for the VBS- 
2 here could be the continuity to the strong VBS-2 phase in 
the broader phase diagram Fig. 2. As is known, the region 
J2 ^ 0.4—2 along the K^iag = axis has strong VBS-2 order 
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FIG. 23: (Color online) Spin and dimer structure factors at a point in 
the VBS-2 phase, Knng = Ji = 1, J2 = 4, measured in the DMRG 
for system size L = 192. The exhibited trial wavefunction is the 
Gutzwiller projection of two equal Fermi seas, {Ni,N2) — (48,48). 
This wavefunction gives decoupled legs expected in the large J2 
limit and does not have any VBS-2 order; however, it reproduces 
the DMRG data quite well, so at this A'ling = Ji cut the system is 
close to the decoupled legs limit even just outside the SBM. 



However, this does not preclude possibility of more phases in 
the model with ring exchanges. Thus, along the way at points 
like A'ring = 0.3, J2 = 1.5 and ifring ~ 0.2, J2 ~ 1.2 we 
also see a dimer feature at q = 7r/2 in addition to a likely 
Bragg peak at q = tt. Since our primary interest is the SBM 
phase, we do not explore the states at large J2 further, loosely 
referring to all of them as VBS-2 in Fig. 2. 
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